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Abstract 

We consider an impurity (AT— level atom) driven by monochromatic light in a host environment which 
is a fermionic thermal reservoir. The external light source is a time-periodic perturbation of the atomic 
Hamiltonian stimulating transitions between two atomic energy levels E\ and Em and thus acts as an 
optical pump. The purpose of the present work is the analysis of the effective atomic dynamics resulting 
from the full microscopic time-evolution of the compound system. We prove, in particular, that the atomic 
dynamics of population relaxes for large times to a quasi-stationary state, that is, a stationary state up 
to small oscillations driven by the external light source. This state turns out to be uniquely determined 
by a balance condition. The latter is related to "generalized Einstein relations" of spontaneous/stimulated 
emission/absorption rates, which are conceptually similar to the phenomenological relations derived by 
Einstein in 1916. As an application we show from quantum mechanical first principles how an inversion of 
population of energy levels of an impurity in a crystal can appear. Our results are based on the spectral 
analysis of the generator of the evolution semigroup related to a non-autonomous Cauchy problem effectively 
describing the atomic dynamics. 

1. Introduction 

In the present paper and in a companion one pQ we study the dynamics of an impurity (an atom) in a crystal, 
or host environment, interacting with an external monochromatic light source serving as an optical pump. 
The host environment corresponds to free electrons in thermal equilibrium within the crystal whereas the 
atom is described by a iV-level system, the pure states of which are the unit vectors of the finite dimensional 
Hilbert space C d (d > N > 1). The external monochromatic light source is a time-periodic classical field 
stimulating transitions between the two energy levels E\ and En- The microscopic non-autonomous dynamics 
of the full system, that is, impurity, host environment and external light source, is then described through a 
two-parameter group of automorphisms generated by a time-dependent symmetric derivation acting on the 
C* -algebra of observables of the compound system. This dynamics is generally non-unitary when restricted 
to the atomic subalgebra. The restricted dynamics of the impurity shows a dissipative behavior, provided the 
coupling to the host environment (thermal reservoir) is effective. In the second part [T] of the present work 
we prove that the dynamics of the impurity is properly described - up to small corrections for a moderate 
reservoir-atom-pump interaction - by some effective non-autonomous time-evolution involving atomic degrees 
of freedom only. See Theorem |3.3| 

The detailed study of this effective dynamics is the main goal of this paper. More precisely, we are interested 
in an effective description of the time-evolution of populations of the atomic energy levels, i.e., an effective 
atomic block-diagonal dynamics. To this end, we use evolution semigroup techniques transforming the non- 
autonomous Cauchy problem into an autonomous one on a suitable Banach space. The same method is also 
used in [T] to analyze the full microscopic dynamics. In fact, many key arguments of the present analysis recur 
in pQ , albeit technically more involved. We remark that evolution semigroups to study the long-time behaviour 
of quantum systems have also been used by Abou-Salem and Frohlich [2]. 
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Observe that similar models without external light source have been extensively analyzed, see for instance 
p2 [H El E] and references therein. One of the most important questions considered previously (see, e.g., [H03 
H]) are existence and asymptotic stability of stationary states, especially asymptotic stability of thermal 
equilibrium states. The latter refers to the so-called "return to equilibrium" which typically occurs in models 
involving one thermal reservoir at fixed temperature weakly coupled to an isolated atom. Note also that, as 
soon as there are several thermal reservoirs at distinct temperatures, the system does not possess a thermal 
equilibrium state, but rather a "non-equilibrium stationary state" (NESS) [HI [TO]. In our setting we certainly 
cannot not expect a thermal equilibrium state to exist due to the optical pump. 

However, it turns out that the atomic block-diagonal dynamics attains for large times a block-diagonal, quasi- 
stationary state g, that is, a stationary state up to small oscillations driven by the optical pump. Observables 
describing quantum coherences between atomic energy levels are not considered. We show that the structure of 
this quasi-stationary state is uniquely characterized by a rather intuitive balance condition involving population 
densities only (Theorem 4.7 1. This balance condition reads (A + B)g = 0. The operator A does not depend on 
the pump and can be interpreted as a matrix of spontaneous transition rates, whereas B is the matrix related 
to the stimulated transition rates, which are induced by the optical pump. We show that A and B satisfy some 
relations which are conceptually related to the well-known Einstein relations [llj . 

Einstein considered an atom interacting with a black body (i.e., broad band) radiation field and derived by 
phenomenological arguments that the stimulated transition rates B between atomic energy levels are propor- 
tional to the radiation density as well as to the spontaneous transition rates A. In contrast, we consider the 
situation of a narrow band (i.e., monochromatic field) driving an atom in interaction with a thermal reservoir. 
In this case, we find that the effective stimulated processes are proportional to the intensity of the radiation in 
this model and that, apart from this simple dependency on the pump intensity, the natural decay rates of the 
atom (i.e., the dynamical properties of the atom as the pump is turned off) uniquely determines the stimulated 
transition rates B. The relevant dynamical parameters turn out to be the decoherence times of the atom, which 
are purely quantum mechanical objects. In a recent work |12j . Berman, Merkli and Sigal show for a similar 
class of models (without pump term, however) that quantum coherences of the atom decay exponentially fast 
for well-defined time scales td cc of decoherence. We show that the (pump independent) parameters td oc and 
the intensity of the pump uniquely define effective stimulated rates. 

As a result of our analysis we find that the time-evolution of the population density of the atomic energy levels, 
i.e., the block-diagonal dynamics, differs from the phenomenological time-evolution usually used in the physics 
literature (cf. [3]). The block-diagonal evolution is usually described by the so-called Pauli master equation, 
which is a phenomenological first order autonomous linear differential equation. Our analysis shows that the 
effective block-diagonal evolution derived from the microscopic model is well described by a integro-differential 



equation (Theorems 4.6 and 4.8). Nevertheless, the long-time behavior of the solutions of our equation is in 
generic situations the same as the one predicted by the Pauli master equation. 

The model we study here is a basic model for pumping schemes of doped crystals as explained in [TOJ. One 
expects for large times a steady emission of photons with fixed frequency E^" — En> , where E\ < ■ ■ ■ < En> < 
■■■ < Sjv" < • • ■ < En are the (ordered) eigenvalues of the atomic Hamiltonian, i.e., the atomic energy levels. 
Such a steady emission should result from a stronger occupation of the energy atomic level En" as compared 
to Ejsf < En». This effect, called inversion of population, is a central mechanism to obtain lasing materials (cf. 
|14j ) and we show that the quasi-stationary state g characterized by the balance condition may exhibit such an 
inversion of population. 

In another recent work [15] . Bach, Merkli, Pedra and Sigal investigate the possibility to control the decoherence 
time(s) tdcc using a different external light source, which does not act as an optical pump but imposes oscillations 
of the atomic energy levels. In certain situations control of decoherence could thus enhance the inversion of the 
population. Consequently, measuring the threshold of the pump intensity needed for inversion of population 
could yield a simple experimental test for control of decoherence. 

To our knowledge there is only one framework in which some aspect of laser phenomenology has been rigorously 
analyzed from first principles, namely for some versions of the Dicke model [TO] , see [TOJ [TOJ [TOJ [20] ■ In Q3 HH] > 
Hepp and Lieb study an interacting conservative system consisting of a reservoir (a radiation field), a finite 
number N of two-level atoms, and finitely many quantized oscillators as "radiation-modes" . In [TOJ [TO] the 
reservoir is absent in their study. In [TOJ [TO] it is shown that this system undergoes a phase transition in the 
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limit N — > oo with the appearance of a coherent radiation driven by the reservoir, whereas (TSJ [19] establish a 
transition from a normal to a supcrradiant phase as N — > oo. More than 20 years later, Alii and Sewell use in 
[213] a dissipative (i.e., non-conservative) version of the Dicke model to get similar results as N — > oo. For more 
details, we recommend for instance [STJ Chap. 11]. In a more recent paper [32] the dissipative model of [2D] 
is proven to be the Markovian approximation of the model considered in |17j . On the other hand, solid state 
lasers are usually constructed with weakly doped crystals [13] . Consequently, an infinite number of impurities is 
not fundamental for the inversion of population. Moreover, the phenomenology of lasers as described in physics 
textbooks is based on three- or four-level atoms [13], but Dicke-type models are based on two-level atoms 
which cannot explain the inversion of population at finite number of impurities. The assumption of impurities 
with three or four levels is also very realistic from the experimental point of view, as most of lasing materials 
used in the praxis are of this type [13] . We focus on a comprehension of the pumping scheme in lasers which is 
(i) coherent with physics textbooks and experimental facts and (ii) uses first principles of quantum mechanics 
only. An important open problem remains however to find a realistic description of a cavity. A recent work in 
this direction is due to Bruneau and Pillet 1231. 



To resume, we focus here on the derivation and analysis of the structure of the quasi-stationary state of 
the corresponding time-dependent master equation for the impurity, whereas its link to the full microscopic 
dynamics is established in [T]. 

The paper is organized as follows. In Section [2] we introduce the microscopic model. Then, in Section [3] 
we define the effective master equation and specify its relation - proven in pQ - to the microscopic model. 
In Section |4.2[ we introduce the evolution semigroup of the non-autonomous effective master equation on a 
suitable enlarged Hilbert space. This leads to a study of an autonomous problem. The spectral analysis of 
the generator of such evolution semigroup performed in Sections |4.3f[4.5| yields a pre-master equation and a 
balance condition, which characterizes the quasi-stationary state for large times in terms of generalized Einstein 
relations as explained in Section [5] Numerical simulations show in Section 4.1 an inversion of population with a 



dynamics far from the one described by the usual Pauli equations. Finally, Section [6] is an appendix consisting 
of a collection of some notions and results on completely positive (CP) semigroups used in our proofs. Moreover, 
it provides heuristics of the proof of the result of [I], some technical proofs as well as a detailed discussion of 
properties of the operator B appearing in the balance condition, which is in general not a generator of a CP 
semigroup. 

Notation 1.1 (Generic constants) 

To simplify our notation, we denote by C, c, c' any generic positive and finite constant. Note that these constants 
do not need be the same from one statement to another. 



2. The host environment-impurity-light source microscopic model 

Laser devices are based on several physical processes. One of them consists of pumping the electron densities 
of atomic energy levels to obtain a so-called inversion of population. It is achieved in many cases by using an 
external source of light. This procedure is called here optical pumping. Note however that pumping can also 
be implemented by means of chemical reactions, electric currents, and other methods. A further step is to use 
this inversion of population to obtain optical amplification through stimulated emission of photons. These two 
processes implement a gain medium, which is then put into a resonant optical cavity to obtain a laser, i.e., 
light amplification by stimulated emission of radiation. Here, we are only interested in the optical pumping 
performed in the gain medium. 

Since the invention of lasers in 1958 by the physicists Schawlow, Townes and Bassov, a significant amount of 
gain mediums have been found. This includes semiconductors, liquids with dyes, gases such as carbon dioxide, 
argon, etc., and solids such as doped crystals and glasses. A large part of gain mediums are based on a host 
material, usually a solid or a liquid, containing impurities or dyes, which can be pumped. A typical example is 
given by the well-known Nd:YAG laser which uses a crystal doped with neodymium. Inversion of population 
is obtained in this family of laser devices by optical pumping. Such a gain medium is typically what we have in 
mind here. 
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An appropriately chosen impurity is crucial to get inversion of population by pumping and hence a positive 
optical gain. In particular, the pumping step should involve at least three or four atomic energy levels, as 
explained in any physics textbook on lasers. For instance, in the case of a monochromatic pump, neodymium 
impurities in Nd:YAG lasers are well described as four level systems. 

Surprisingly, also dissipative processes are important in order to get an efficient optical pumping, in general. 
For instance, we could analyze the model for impurities in crystals which will be defined below with A = 
(meaning that the impurity does not interact with the surrounding environment) and r\ ^ (meaning that the 
optical pump is turned on). For such a choice of parameters, the dynamics can be explicitly computed, showing 
that the atomic populations undergo Rabi oscillations and one cannot obtain a (quasi-) stable inversion of 
population and not even a positive inversion in time average. In solid-state lasers, the interaction with the host 
material, i.e., the crystal or the glass, provides a dissipative component to the (effective) atomic dynamics. The 
same should be found in a gain medium constituted of a liquid with a dye to be optically pumped. 

The dissipative mechanism can be quite different from one gain medium to another. It can be due to many 
different types of interactions such as interactions with phonons, electrons in a crystal, etc. Nevertheless, its 
only important property for the pumping process is to damp Rabi oscillations of pumped atoms and to modify, 
in interplay with the optical pump, transition rates between atomic energy levels in order to get a (quasi-) 
stable inversion of population. 

Therefore, we can conclude three things about the optical pumping, i.e., the stimulated process to get an 
inversion of population: 

(a) To explain optical pumping, it suffices to use a generic dissipative mechanism. To this end, we consider 
here an interaction between a N— level atom (serving as a small system in a host environment) and a 
fcrmionic thermal (macroscopic) reservoir whose mathematical setting is standard, see, e.g., [531 123- This 
choice excludes interactions between impurities and bosonic particles like phonons. This could also be 
implemented, but for simplicity we refrain from considering it as it qualitatively leads to similar results. 
Note, however, that more physical models for host environments could be useful to get more precise 
quantitative results. 

(b) To obtain an efficient optical pumping, the strengths r), A > of respectively the optical p ump and the 



interaction with the host environment have to be of the same order, i.e., rj ~ A, see Section 2.7 Indeed, 
if 77 >> A, then Rabi oscillations with frequency of order fj are the dominating processes governing the 
dynamics of populations. This oscillations are then generally damped in a time-scale of order A not 
depending much on the intensity of the pump. In the opposite situation, i.e., when 77 << A, no inversion 
of population can appear as the system relaxes in this case to a state near the ground state (or a thermal 
equilibrium state) of the atom and all the energy provided by the (weak) optical pump is lost into the 
(very large) host environment. 

(c) Finally, the inversion of population described in physics textbooks makes sense if the "TV-level atom" 
picture stays valid. In other words, the interaction with the host environment must be a small perturbation 
of the Hamiltonian representing the A-level atom. By (b), the optical pump should also be seen as a 
perturbation of the N— level atom. In view of this observation it is appropriate to use (Kato's) perturbation 
theory 26J of (discrete) eigenvalues of the A-level atom. 

The host environment discussed in (a), i.e., the fermionic thermal reservoir, is described in details in Section 
2.1| We then define the impurity and the external monochromatic light source in Sections |2.2[^2.3| respectively. 



The coupled full system is set up in Sections 2.4 2.7 



2.1 The host environment as a fermionic thermal reservoir 

Let f)i := L 2 (M 3 , C) be the separable Hilbert space representing the one-particle space of the host environment 
(reservoir). The one-particle Hamiltonian hi is then defined by using a dispersion relation. Indeed, we choose 
some measurable, rotationally invariant function E : M 3 —> K, i.e., E(p) = E(|p|), and define the multiplication 
operator hi = ft-i(E) by f(p) 1— > E(p)/(p) on (ji. Physically, E(p) represents the energy of one particle with 
momentum p within the reservoir. 
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In the case where the host material of the gain medium is a crystal, observe that L 2 spaces on Brillouin 
zones as one-particle spaces are more realistic than i 2 (]R 3 ,C). As explained above (cf. Observation (a)), the 
results obtained would qualitatively be the same. For instance, the rotation symmetry assumed above for the 
dispersion relation E is not essential to the analysis performed below. However, to satisfy the assumptions of 



Theorem 3.3 we require that E behaves (at least near p = 0) like \p\ up to some diffeomorphism. For instance, 
the dispersion relation E(p) = \p\ 2 is allowed as \p\ and \p\ 2 are clearly the same function up a diffeomorphism 
ofM 3 \{0}. 

The field algebra Viz of the reservoir is the CAR C*-algebra generated by the annihilation and creation 
operators a(/), a + (/) := a(/)*, / <E F)i, acting on the antisymmetric Fock space .F_(hi) and satisfying the 
canonical anti-commutation relations (CAR): 

a{h)a+{f 2 ) + a + (f 2 )a(h) = (/ lf / 2 ) , a(h)a(f 2 ) + a(f 2 )a(h) = 

for any fi,f% 6 rji- (•, •) denotes the scalar product in Ijj. 

The (unperturbed) dynamics of the reservoir is given by the family {if } t6 K of Bogoliubov automorphisms 
on the algebra Viz uniquely defined by the condition: 

if (a(/)) = a(e 4t/l1 /), / € f)i, teR. (2.1) 

Physically, this means that the fermionic particles of the reservoir do not interact with each other, i.e., they form 
an ideal Fermi gas. The group t u :— {if } tS R of automorphisms is strongly continuous and hence, (VtzjT 71 ) is 
a C*-dynamical system. We denote its generator by 8-r- 

Note that generators of C*-dynamical systems are symmetric derivations. This means that the domain 
Dom(<5-ft) of the generator S-r is a dense sub-*-algebra of Viz an d, for all A,Be Dom(S-ji), 

Stz(A)* = S n (A*), Sn(AB) = 5 n (A)B + A5 n {B) . 

Thermal equilibrium states of the reservoir are defined through the bounded positive operators 
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acting on f)i for all inverse temperatures [3 e (0,oo). Indeed, the so-called symbol d-jz uniquely defines a 
(faithful) quasi-free state 

utz ■■= Lo d ^ (2.2) 
on the fermion algebra Vtz by the conditions Wd^i^-n) = 1 and 

cod-*. ( a+ (fi) ■ ■ ■ a + (/m)a(ffi) • ■ ■ a(g n )) = 5 m . n det ([(gj, d n fk}]j,k) 

for all {/j}" =1 , {<7j}™ =1 C f)i- The positive, normalized linear functional wjj is the unique /3-KMS state of the 
C*-dynamical system (V-r,, t k ) and is called the thermal equilibrium state of the reservoir at inverse temperature 
/3 6(0, oo). 

This definition of thermal states is rather abstract, but it can physically be motivated as follows: Confining 
the particles within a box of side length L corresponds to the replacement of the momentum space M 3 by ^Z 3 , 
i.e., L 2 (K 3 ) by ^ 2 ( 2 ^ L Z 3 ). In particular, the spectrum of hi (and hence of its fermionic second quantization 
dT^(hi)) becomes purely discrete. Additionally, the operators e _ ^ dr -^ 1 ^ are in this case trace-class for all side 
lengths L. Hence, we can define Gibbs states 

a W n ._ Tr(. e -^-(M) 
»k U • Tr ( e -/3dr_(M) ' 

which have the thermal equilibrium state lu-jz as unique weak-* limiiQfor L — > oo. This follows, for instance, 
from the results of [551 Chapters 5.2 and 5.3] on KMS states. 

1 This refers to the weak* -topology on the locally convex real space V^, the dual space of the separable Banach space V-n- 
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2.2 The impurity as a N level atom 

The impurity (atom) is modeled by a finite quantum system, i.e., its observables are the self-adjoint elements 
of the finite dimensional C*-algebra B(C d ) of all linear operators on <C d for d G N (d > 2). 

In the sequel it is convenient to define left and right multiplication operators on B(C d ): For any A G B(C d ) 
we define the linear operators and ^4 acting on B(C d ) by 

B i-> AB := AB and B i-> ^45 := . (2.3) 

The Hamiltonian of the atom is an arbitrary observable H at = i?* t G B(C d ) representing its total energy. 
We denote its eigenvalues and corresponding eigenspaces by Ek G K and %k C C d for fc G {1, . . . , A 7 "} (AT > 2), 
respectively. Ek is chosen such that Ej < Ek whenever j < fc. In other words, Ek is the energy of the fcth atomic 
level and vectors of T-Lk describe the sub-band structure of the corresponding energy level. The dimension rik 
of the eigenspace Hk is the degeneracy of the fcth atomic level. 

As usual, the Hamiltonian H at defines a free atomic dynamics, i.e., a continuous one-parameter group of 
automorphisms r at :— {T at } tS R of the C*-algebra £>(C d ) defined by 

rf{A) := e ttH ^Ae~ itH ^ , A G B{C d ) , (2.4) 

for all t G M. 

Thermal equilibrium states of the free atom are Gibbs states g a t given by the density matrix 

g— 

P * := Tr c , ( e -^.t) ( 2 ' 5 ) 

for any inverse temperature (3 G (0, oo). 

In presence of the optical pump and the host environment (the thermal reservoir), the state of the atom is 
generally far from the Gibbs state g at . We thus consider arbitrary atomic states w at . For any state w at on 
B(C d ), there is a unique trace-one positive operator p at on C d , the so-called density matrix of cj at , such that 

w at(^4) = Tr C d (p at A) , AeB(C d ) . 

Note that any state w a t on B(C d ) can be represented as a vector state via its GNS representation (J5at ; ""at, ^at), 
see, e.g., [SU Theorem 2.3.16]. If w at is faithful then (io a t, 7r at , O a t) is explicitly given as follows. The Hilbcrt 
space .f) a t corresponds to the linear space B(C d ) endowed with the Hilbert-Schmidt scalar product 

(A,B) at := Tv cd (A*B) , A, Be B(C d ) . (2.6) 



The representation 7r at is the left multiplication explained above in (2.3|, i.e., 

7r at (A) = 4. , £(C d ) . 

The cyclic vector of the GNS representation of w a t is defined by using the density matrix p at G £?(C d ) of Lu at as 

n at := pll 2 G j5. t . (2.7) 
Using the cyclicity of the trace we obtain that 

Wat(A) = <n at> 4fl at > at , ^G6(C d ). 

This GNS representation (i^ a t, 7r at , O at ) is known in the literature as the standard representation of the state 
w at - See [13 Section 5.4]. 

The dynamics given by the continuous one-parameter group r at of automorphisms of the C*-algebra B(C d ) 



defined by (2.4) can be represented in the Schrodinger picture of Quantum Mechanics through the so-called 
(standard) Liouvillean operator 

L a t '■= H a ^ — H at = [H at , ■ ] = L* t (2.8) 
acting on the Hilbert space ^ at . Indeed, it is easy to check that: 
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Lemma 2.1 (Schrodinger picture of r at ) 
For all t £ R, 

CJ at (rf (A)) = (0 at (t) , 7T at (A) fi at (t)) at , A G £(C rf ) , 

w/iere fi at (i) := e- itL »*fi at . 

Finally, for fc G {1, . . . , N}, note that the population of the fcth atomic level in the state w at is defined by the 
expectation 

p fc (p at ) := « at (1 [^at = E k ]) = Tr c , (1 [H at = E k ] p a% ) > , (2.9) 

where 1 [H at = iSfc] G £>(C d ) is the orthogonal projection onto the eigenspace Hk- If w a t = 0at is the Gibbs state 
of the atom and the atomic energy levels are non-degenerated or the inverse temperature /3 G (0, oo) is large 
enough then 

Vj, k € {1, ...,N} , j < k : P] (p g ) > Pk (p B ) ■ 

In contrast, we say that a state w or a density matrix p shows inversion of population if there are j, k G {1, . . . , N} 
such that j < k, i.e., Ej < E k , and Pj(p) < Pk{p)- In other words, inversion of population requires a higher 
energy level more populated than a lower one. Of course, this phenomenon can, in general, only appear in a 
state out of equilibrium and one usually uses external light sources to artificially pump electrons from a low 
energy level of the atom to a higher one. 



2.3 The external monochromatic light source as a classical optical pump 

The optical pump, i.e., the monochromatic photon field interacting with the atom, is described by the following 
time-periodic perturbation of the atomic Hamiltonian H at : 

r]cos(mt)H p , tjj :— En — E\ > , tel. (2.10) 

Recall that E\ < ■ ■ ■ < E N denote the N eigenvalues of H at = H* t G B(C d ). Here, 

H p := hp + h* p G B{C d ) (2.11) 

for some h p G B(C d ) satisfying 

ker (hp) 1 - C U x := ran (1 [H at = E{\) , (2.12) 

ran(/i p ) C H N := ran(l [H at = E N }) . (2.13) 

In other words, the optical pump produces only transitions between the lowest and the highest atomic levels 1 
and N, as described in standard textbooks on the physics of lasers. 

From the physical point of view, the time-dependent optical pump may be regarded as a partial classical limit of 
a closed (autonomous) physical system involving a quantized pump modeled by a quantum harmonic oscillator. 
The corresponding initial state for this quantized pump should be chosen as being a coherent state. See, e.g., 



Remark 2.1 (Non monochromatic light sources as classical optical pumps) 

Results of this paper can easily be extended to non-monochromatic light 
sources as classical optical pumps. This case corresponds here to replace the cosine in \2.10 ) by 



some 



time-periodic and continuous function. However, in order to keep technical aspects as simple as possible, we 
refrain from considering this more general case. 



2.4 The uncoupled reservoir atom system 

Define the C*-algebra V := B(C d ) ® V-r. As both C*-algebras B(C d ) and Vn are already realized as algebras 
of bounded operators on Hilbert spaces and since B(C d ) is finite dimensional, we do not have to specify the 
meaning of the tensor product. Observables of the reservoir-atom system are self-adjoint elements of V. Its 
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free dynamics is described by the strongly continuous one-parameter group r := {r t }teR of automorphisms of 
V defined by 

Tt :=rf ®rf , tel. (2.14) 

This tensor product is well-defined and unique because the atomic algebra B(C d ) is finite dimensional. The 
generator of the free dynamics defined by r is a symmetric derivation, denoted by 6, which acts on a dense 
sub-*-algebra Dom(£) of V. 

Let w a t be any initial (not necessarily Gibbs) state of the atom and define the initial state of the atom-reservoir 
system by 

w := w at ® w-r, . (2-15) 

Again, the latter is well-defined and unique, by finite dimensionality of B(C d ). If w at = g a t is the Gibbs state 
then w is clearly a (j3, r)-KMS state. Observe also that g at is a faithful state and we assume without loss of 
generality that w at is also a faithful state. Indeed, the set of faithful states is dense in the set of all states of 
the atom. Since the quasi-free state wr of the reservoir is also faithful, this property carries over to the initial 
state wq of the composite system. 



Remark 2.2 (Coupled initial state of the atom reservoir system) 

Considering impurities interacting with the host environment long before the pump is turned on, the initial state 
of the atom-reservoir system should, in principle, not be a product state as in (2.15). Instead, it should be a 
thermal equilibrium state of the coupled atom-reservoir system. However, in rather generic situations it can 
be shown that at small atom-reservoir couplings this thermal state is near the product state wo with w at = g a t 
and the results would be the same up to sub-leading corrections. As the KMS states of the model considered 
here are unique (KMS states of bounded perturbations of a free fermion gas are unique), the latter follows from 
standard results on the stability of KMS states, see for instance 1251 Section 5.4-. 1]. Indeed, we can even treat, 
by the same methods, any initial state of the composite system as soon as its relative entropy with respect to the 
product state (2.15) is finite. This feature is verified for the equilibrium thermal state of the composite system 
at weak coupling. In order to keep technical aspects as simple as possible, we will not consider this case. 



2.5 The atom reservoir interaction 

The interaction between the atom and the fermionic thermal reservoir involves the so-called fermionic field 
operators defined, for all / £ fji, by 

*(/) := + °(/)) = e B (-M<h)) ■ 

Choose now a finite collection {Q e }^ =1 C B(C d ) of self-adjoint operators and an orthonormal (finite) system 
{ftYeLi C Then, the atom-reservoir interaction is implemented by the bounded symmetric derivation 

m 

e=i 



Note that the orthonormality of the family {fe}JLi does not inflict loss of generality. Indeed, for an arbitrary 
finite set {Q e }^L 1 C B(C d ) of self-adjoint operators and (possibly not orthonormal) family {fe}"Li ^ 
there are m £ N, a finite collection {Q e }™ =1 C B(C d ) of self-adjoint operators and an orthonormal system 
{fe}"Li c &i such that 

m m 

i=i e=i 

As we seek to maintain mathematical rigor while keeping technical aspects as simple as possible, we assume 

to 



some technically useful conditions on the family {fe} i=1 C f)i guaranteeing the assumptions of Theorem 3.3 
be satisfied. Note that these conditions will by no means restrict the range of our analysis to rather physically 
meaningless submodels. First, {fiYi^i C t)i is taken as a family of rotationally invariant functions, i.e., fi{p) = 
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fiflpl) for all pgK 3 and I G {1, . . . , m} with : Rj — > C. The rotational invariance of for all £ G {1, . . . , m} 
is only assumed for technical simplicity as more general choices of such functions would yield similar results, at 
least qualitatively. The second, rather technical assumption is an analyticity condition which is only required 
to prove Theorem |3.3| It is given here for completeness, but it can clearly be omitted as no proof of this paper 
uses it. This analyticity condition reads as follows: For all £ G {1, . . . , m}, the complex valued functions gi and 
gf respectively defined by 

_ i/9 r k(x) , x>o, 

VzgR: ge(x) := |ar| (l + e^ x ) 1 I (2.16) 

{ h{-x) , x < , 

and gf(x) := ige(-x) on R have an analytic continuation to the strip R+i(— C, C), and satisfy 



sup \ / (\ge(x + \gf(x + i-&)\) 2 dx } < oo 

#e(-C,C) Ui J 

for all I G {1, . . . , m}. For instance, to satisfy these conditions one may choose for any I G {1, . . . , m}, the 
function as linear combinations of terms of the form |a;| 2 ' c_1 exp(— cx 2 ) with k G No- 

2.6 Dynamics of the coupled atom reservoir pump system 

The full dynamics of the system involves the classical pump described in Section |2.3[ which is implemented 



as a periodic perturbation of the dynamics of the atom-reservoir system described in Section 2.4 Indeed, let 

6at, P := i[H p ® l n , ■} 

and rj G R. The coupled atom-reservoir-pump dynamics is then generated by the time-dependent symmetric 
derivation 

S[ x <ri : = 8 + 7]cos(wt)S attP + XS^n, t€l. (2.17) 
Here, A, 77 G R are the atom-reservoir and atom-pump coupling constants, respectively. 

Observe that S 1 ^'^ acts on a dense sub-*-algebra Dom(<5t) = Dom(J) of V which does not depend on t G R. 
Indeed, 5 a t,K and <5 at , p are bounded symmetric derivations and (J^' 1 ^ is the generator of a strongly continuous 
one-parameter group of automorphisms of V. As the map 

1 1 — ^ S^ v) - 5 ( X ' V) 

is norm-continuous, 5\ X ^ generates a strongly continuous two-parameter 
family Wt^^tys of automorphisms of V corresponding to the non- 
autonomous dynamics of the (coupled) atom-reservoir-pump system. The operator can even be 
explicitly constructed as a Dyson series because the operator r/ cos(mt)5&t,p is bounded and cos is a smooth 
function. 

The time-evolution of the state of the full system is then given by 

0J t := U O Ti^ V ' = (w a t ® OJ-n) O Tip 71 ', t G Rq" . 

The reduction of this state onto the atomic degrees of freedom yields a time-dependent atomic state defined by 

w. t (t)OA):=wt(.A®lw), AGS(C d ), (2.18) 

for all < G Rq~. 

2.7 Moderate optical pump and atom reservoir interaction 

As explained at the beginning of Section [2] (cf. Observation (c)), we are interested in the regime where 
|A|, \r]\ << 1. In other words, we take the atom-reservoir and atom-pump interactions as being small, but non- 
vanishing, perturbations of the free dynamics generated by the symmetric derivation S. Moreover, we assume 
that the pump is moderate with respect to the atom-reservoir interaction in the following sense: 
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Assumption 1 (Moderate optical pump) 

For any Asl and some fixed constant C £ (0, oo): \rj\ < CX . 

Actually, in all our proofs, it would suffice to impose \n\ < C\X\ for some sufficiently small constant C £ (0, oo). 



However, as it will be shown below (cf. Remark 4.3), the contribution of the pump to the final atomic state is 
of order rj 2 /X whereas the contribution of the interaction with the reservoir is of order one (in the parameter 
ry 2 /A 4 ). Thus, imposing \rj\ ~ A 2 means physically that both the pump and the reservoir contribute in an 
essential way to the final state of the atom. The (relative) strengths rj, A > of respectively the optical pump 
and the interaction with the host environment turn out to be equal to fj :— n 2 /X 2 and A := A 2 . In particular, 
Assumption [l] means that fj ~ A. See Observation (b) above. Consequently, we say in this context that the 
pump is weak whenever \rj\ « A 2 , i.e., fj << A. 

In the opposite situation when for small \rj\, |A|, (77, A) £ K 2 , one has \n\ » A 2 , i.e., fj >> A, Rabi oscillations 
are the d omin ant processes driving the populations of the atomic energy levels. Indeed, a similar version of 
Corollary 4.5 is still valid in the strong pump regime \n\ » X 2 if \rj\ « \X\. Using this result one can show that, 
in general, Rabi oscillations dominate the dynamics of populations at time-scales comparable to and that 
there is a global relaxation of populations at time-scale X 2 not depending much on the pump intensity. By 
this last property we could also call the regime \n\ >> X saturated pumping. However, we are rather interested 
in the situation where pump and natural relaxations of the atom compete in a non trivial way with each other 
and the evolution of the full system is well described by some relaxing dynamics driving the atom to a state 
with a persisting inversion of population. See again the discussions at the beginning of Section [2j in particular 
Observation (b). 



3. The effective atomic master equation 
The aim of this paper is to analyze the atomic dynamics resulting from the r estriction on B(C d ) of the 



full dynamics generated by the symmetric derivation 6 



(A,»j) 



family {u at (t)} t M + of states defined by (2.18) or, equivalently, to the family {p at (t)} 



See Equation (2.17). This corresponds to the 



> le 



E + of density matrices. 



More precisely, we are interested in the timc-bchavior of observables related to atomic levels only, and not to 
correlations between different levels. This amounts to study the orthogonal projection (p at (£)) of the atomic 
density matrix p at (t) on the subspace 



S = 2)(Jf at ) :=B(W, 



(3.1) 



of block-diagonal matrices. In other words, we analyze the density matrix 

N 



Pv (Pat (*)) = E 1 = E k ] Pat (t) 1 [H at = E k ] 



(3.2) 



k=l 



for any teKj. 

As proven in pQ, the density matrix p &t (t) can be well approximated, uniformly in time, on the subspace 2D 
by the solution of an effective non-autonomous initial value problem in B(C d ), the effective atomic master (or 
Kossakowski-Lindblad) equation described in Section 3.3 Th e corresponding time evolution is generated by 



the time-dependent Lindbladian defined in Section 



and its derivation from the microscopic model is 
discussed in 6.1 In Section 3.2 we impose a condition, introduced in [301 Theorem 2] in the context of relaxing 



3.1 



CP semigroups, on the dissipative part of this generator which is a non-commutative version of the irreducibility 
of classical Markov chains and ensures the existence of the (quasi-) stationary state for populations at large 
times. 



3.1 The atomic time dependent Lindbladian 

Lindblad operators (or Lindbladians) are standard objects used to describe (generally dissipative) dynamics 
compatible with some phenomenologically reasonable prescriptions like complete positivity. As explained in 
[5J , the typical example of application of such operators is related to the time-evolution of composite systems 
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made of a small (open) quantum system interacting with a macroscopic one (a reservoir). In this context, the 
dynamics restricted to the small system is in many situations well-described by a semigroup {e t£ } t>0 which 



is trace-preserving and completely positive, see Section 6.3 Generators £ of such completely positive (CP) 
semigroups are called Lindblad operators or Lindbladians. 

A first rigorous result in this direction is due to Davies [31] |32j [33] in certain scaling limit, the so-called 
weak coupling limit for similar interacting systems as ours with n = 0, i.e., without the optical pump or any 
other time dependent term. See also [S7J and references therein. For more details on CP semigroups, we also 
recommend Section |6.3[ Note however that, in contrast to Davies' approach, we never take the limit A — > 0. 
It suffices in our analysis to have a sufficiently small coupling constant |A| 7^ 0. Observe also that we have to 
control at the same time the (non-trivial) influence of two microscopic couplings of different orders (\r]\ ~ A , 
by the moderate pump assumption) and the precise meaning of a weak-coupling limit is not clear from the 
beginning. The dynamics we deal with is moreover non-autonomous. 

The physical system considered here yields a non-autonomous effective time-evolution with a time-dependent, 
generally non anti-self-adjoint generator £j A:I? ' € B(Sat) [J given, for any t £ R, by: 

£ { t X - v) {p) :=£ at (p)+77Cos(tui)£ p (p) + A 2 £ K (p), p £ Sat . (3.3) 

See Section |6.1| for an heuristic derivation of this time-dependent Lindbladian from the microscopic model. 

The first term defining £,[ X ' r '^ is the Lindbladian of the free atomic dynamics. It is the anti-self-adjoint 
operator 

£ at (p) := -iL at (p) = -i[H at ,p] = -2* at (p), p e Sat . (3.4) 



See Lemma 2.1 The second term of (3.3 1 encodes the influence of the optical pump. It is defined via the 
Lindbladian 

£ p (p) := -i[H p , p] = -£* (p), p e Sat • (3.5) 

The third term includes a dissipative part 2d & #(Sat) corresponding to the non-unitary character of the 
effective dynamics and so, 2,-r € £>(Sat) is not anti-self-adjoint. In fact, the Lindbladian £r, is related to 
the second-order perturbation term coming from the atom-reservoir (electron field-impurity) interaction in a 
suitable representation and equals 

£k(p) ■= -i[H Lamh , p] + £d(p), p £ Sat • (3.6) 

In order to define the so-called atomic Lamb shift -ffLambj which encodes the shift of the atomic energy levels 
due to the presence of the reservoir, and the effective atomic dissipation £,d some preliminary definitions are 
necessary: We denote the spectrum of any operator A by a (A), its positive part by er + (A) := {A) n Rq~, and 
define the sets 

t e := {(j, k) : E 3 - E k = e} C {1, 2, . . . , N] x {1, 2, . . . , N] (3.7) 

for each eigenvalue 

e e a(L at ) - a([H at , ■ ]) = {E s - E k : j, k G {1, 2, . . . , N}} . (3.8) 
Let {Vj k }j,k,i C B(C d ) be the family of operators defined by 

V$ := 1 [H at = E s ] Q e 1 [H at = E k ] (3.9) 

for j, k E {1,2,..., N} and £ G {1,2,..., m}, and let {f\ 0) }f =l be the family of functions K -)• defined by 

/f'W-^=^feWf (3-10) 



at any fixed inverse temperature f3 £ (0, 00) of the fermionic reservoir, see (2.16). 
Remark 3.1 (Self adjointness of the space spanned by {vj e ^}j t k.e) 

Since Qg = Q\ [, one has (Vj k )* — V}f} f or an U h k G {!) 2, ■ ■ ■ , N} and £ £ {1, 2, . . . , to}, and obviously, 

span{Vj^} jt k,e = span{(V ? ( j. ) )*} j!fc/ c Sat ■ 
This fact is important when using Theorem \6.10[ 
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Then, the atomic Lamb shift ifLamb € B(C d ) is the self-adjoint operator defined by 

h^--=-\ E E E42y#**# ( 3 - n ) 



2 

eetr([H. t ,-])\{0} (i,fc)et« fci 



with the real coefficients 

fW(. + (E k -Ej)) 



d% := VV 



f»j — 

a-d • 2 



being the principal part W[f] of functions / = ff' (• + (f? fe — Ej)). 

Meanwhile, the non anti-self-adjoint operator Zd G £>(-f) a t) describing the effective atomic dissipation is 
defined by 

\ E E E«M2. ( 3 - 12 ) 

eeo-([H» t ,-]) (j,fc)et e <=i 
where c$ := nf { /\E k - Ej) and 

The terms V jk pV jk in these definitions correspond to transitions from the fcth to the jth atomic levels, whereas 
the other terms guarantee the Markov property of the dynamics, i.e., the preservation of the trace of the density 
matrix. 

Note that the functions satisfy the equality 

fW(- x ) =e -P*fV0( x ) i XE R, 

because Q e is 

dissipation equals 



(£) (£)* 

whereas V j k = V k ■ because Q e is self-adjoint, by assumption. Using these properties, the effective atomic 



** = \ E E E c S(42+a-Me-^g) • 

e£<r+([Hat,-]) (j,fe)et_ e <=1 

Observe further that [ifLamb>P fl ] = and 

e »tffat£W e -«^at _ e l * e £^ 

for all (j, fc) G t e , e G u + {[H at . •]). This is the standard form of a Lindbladian fulfilling the so-called (quantum) 
detailed balance condition in the sense of Alicki-Frigerio-Gorini-Kossakowski-Verri with respect to the atomic 
Gibbs state jj at . See [34], [35] which is reviewed in [23 Section 4.5]. One important consequence of this fact is 
that the atomic Gibbs state, for all A, t, satisfies 

4 m) (p b ) = o. 



See, for instance, [3] Section III. 4]. Here, the parameter f3 of the density matrix p g (see (2.5)) is chosen to be 
the inverse temperature of the reservoir. 

Remark 3.2 (Lindbladians as generators of Markov CP semigroups) 

The operators if at , ff p , if Lamb G B(C d ) are self-adjoint and {Cj\}j,k,e are non-negative numbers. Thus, at any 
fixed time t G M and for all (A, rf) G M. 2 . the Lindbladians £j and ^£ p + \ 2 £,-r are generators of Markov CP 
semigroups, by Theorem\6. 



3.2 Irreducibility of quantum Markov chains 

In principle, even after having extracted (by some averaging procedure, for instance) the oscillations of 
frequency w caused by the presence of the pump, the family {P® (p at (i))} tgR + of density matrices could have 
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several accumulation points, limits depending on the initial state, or even be oscillating (Rabi oscillations) as 
t — > oo. We would like to avoid this situation and spectral properties of the Lindbladian ^£ p + A £7?, turn out 
to be important in this sense, see Section |6.2| To this end, we require that is a non-degenerated eigenvalue 
of §£ p + \ 2 £ n with some non-trivial real spectral gap, that is, 

max {Re {w} | w € a (~£ p + \ 2 £ n ) \{0}} < -A 2 C < 

with C € (0, 00) being some fixed constant not depending on A and 77. 



This is useful (and very likely even essential) to prove Theorem 3.3 because it yields uniform bounds in time 



and allows the study of the asymptotic dynamics of the atom. The following assumption on the dissipative part 



£d suffices to ensure the spectral properties mentioned above (cf. Lemma 6.3) 
Assumption 2 (Irreducibility of quantum Marko v ch ains) 



The family {V^}j t k,t C B(C d ) of operators defined by (3.9) satisfies 



( U {*$})" 

{(J,k,e) : c$^0} 



with M" being the bicommutant of M C B(C d ). Recall that cf k := nf^\Ek — Ej), see also (3.10) 



The assumption above highlights the role played by dissipative effects of the fermionic reservoir on the atom in 
order to get an appropriate asymptotic evolution of populations of atomic levels. Actually, the existence and 
uniqueness of the final (t — > 00) density matrix projected on the subspace D of block-diagonal matrices follows 
from this hypothesis (cf. Theorem 4.7). See also Observation (a) at the beginning of Section [2j 



Assumption [2] is a non-commutative version of the irreducibility of classical Markov chains. To illustrate this, 
we consider the following example: Assume for simplicity that m = 1 and the degeneracy of the fcth atomic 
level equals = 1 for all k € {1, . . . , N = d}. Let {fk}t=i <- ^ ^ e an orthonormal basis of eigenvectors of 
-Hat with _ff at (/3 fc = Ek^fk- If the self-adjoint operator Q\ of the atom-reservoir interaction is defined by 



h<Pk 



^2<Pj> fee {!,..., d} 



then the family {V^ k } d k=l satisfies Vj^V^l = vj k for all j, k, n and forms an orthonormal basis of f) a t- In the 
orthonormal basis {<Pk}k=v ^jk ^ s ^ rie elementary matrix made of zeros except at the intersection of the jih row 
with the fcth column where its matrix coefficient is 1. We assume the irreducibility of the family {c^ k } d k=1 C Kq" 
of non-negative numbers in the sense that, for all j ^ k, there is a finite sequence (j'i, ki), . . . , (j n ,k n ) such 
that c^ ki , . . . , c^k n 7^ 0, ji — j, k n — k, and fc; = for I G {1, 2, . . . , n — 1}. Physically speaking it means 
that any arbitrary pair of atomic levels is connected by non-vanishing transitions. By using the commutator 
identity 

frv}$] - = V$[A,V%] + 

for all j, n, k and the irreducibility of the family {c^l}^ k=1 one can compute the commutant 

( U {!$})' -CIc, 

{(J,k) : cWjtO} 

from which Assumption [2] follows. This is in perfect analogy to well-known results about uniqueness of invariant 
states of (aperiodic irreducible) discrete Markov chains. See for instance [Ml Chapter 18]. 

Assumption [2] concludes the list of required conditions and from now on, we assume Assumptions [l}|2] to be 
satisfied. 
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3.3 The effective master equation 

We define now the effective atomic master equation on the Hilbert space Sj at as the initial value problem 



Vt > 



P (t) = ^(Pm P(0) = Pat (0) = Pat € i3a 



(3.15) 



Recall that p at is the density matrix of the initial atomic state w a t of the atom and p at (t) is the density matrix 
of the time-dependent state w a t (t) defined by (2.18) for any t g Kq . Even if one imposes the condition p at g S), 
note that p at (t) is generally not block-diagonal, i.e., p at (t) ^ J) for all t > 0. The same is true for the solution 



of (3.15) 



The effective atomic master equation has a unique solution which, by finite dimensionality of f) at , is explicitly 
given by a Dyson series. In particular, this initial value problem defines a two-parameter family denoted by 



{f( A s ' ,, ''} t > s . Since the Lindbladian S,^' 11 ' is continuous and generates a Markov CP semigroup at any fixed time 
(eR (cf. Remark [3~2| , the two-parameter family {f[ X ^} t > s C B(Sj at ) is continuous, completely positive and 
preserves the trace. The positivity and trace preservation imply that this family is uniformly norm bounded^] 



VA, n, w, s, t G K, t > s 



< C 



(3.16) 



for some finite constant C £ (0, oo) not depending on A, n, w, s, and t. When the optical pump is absent, the 
dynamics becomes autonomous and the family {T[ X g°^}t> s corresponds to an one-parameter semigroup denoted 
for simplicity by 

(3.17) 



,(A,0) 



The main interest of the initial value problem (3.15) is that its (unique) solution p(t) accurately approximates 
; small couplings the true density n 
of block-diagonal matrices for all t € 



at small couplings the true density matrix p at (t) of the time-dependent state ui at (t) on the subspace 2) (3.1 1 

, Indeed, we prove in [1] the following assertion: 



Theorem 3.3 (Validity of the effective atomic master equation) 

Assume that p &t € 2). The unique solution {p(t)} t >o of the effective atomic master equation (3.15) and the 
atomic density matrix {p at (t)} t >o satisfy the bound 

\\Pv (pAt) ~ P(t))\\ <CW|A| 

for some constant C ro € (0,oo) depending on w, but not on the initial state w at of the atom and the parameters 
t, A, and rj. 

Sketch of the proof. The proof of Theorem [33] is conceptually similar to what is done in Section|4j but technically 
much more involved: 



Similar to the one pa ram eter semigroup {T a } a >o defined below, we represent the non-autonomous evolu- 
tion {Jj[ g' }t>s ( CI - (6.2)) as an autonomous dynamics {e ae } Q >o on an enlarged Hilbert space fj 0V o 3 Sj 



of periodic i^-valued functions (vectors of Sj are naturally identified with the constant functions in this 
case). The generator Q of the new time-evolution is occasionally referred to as Howland or Floquet op- 
erator, fj is the GNS-space of the initial state loq and {f/ t g'' 7 "'}^ is a suitable representation of the 



full microscopic dynamics t^q through a family of bounded operators on f). See Section 
explicit construction of Sj and {uj- X s ,r] ^}t> s . 



6.1 



1-4 for the 



Then, we perform an analytic deformation Q{9) (more precisely, analytic translation) of the unbounded 
closed operator Q and prove that the dynamics driven by {e Q ^} a >o and {e"^^^}o,>o are the same on 
the atomic subspace ^ a t = £)at ® {^7?.} C i3 e vo- The use of analytic deformations is the reason for the 



analyticity condition stated at the end of Section 2.5 which is also related to (6.4) 



2 This can be seen by using a decomposition of any A S flat in imaginary and real parts, each of them being also decomposed 
in positive and negative parts. Use then the trace— norm, which is equivalent to the norm on fl a t = B(C d ). 
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• In contrast to Q, whose eigenvalues are all embedded in the continuous spectrum, G(Q) has discrete 
spectrum. It turns out that the discrete eigenspace of Q{6) is the relevant one for the atomic dynamics. 
We then analyze the discrete spectrum and eigenspace of Q{9) through Kato's perturbation theory for 
closed operators. If Assumption [T] holds to the leading order in A and 77, that is, second order in A and first 
order in 77, the operator Q{9) is - up to purely imaginary constants - unitarily equivalent to ^£ p + A 2 £^ 
in finite dimensional invariant subspaces spanning the whole discrete subspace of G(0). 

• We are then in position, by using the inverse Laplace transform for Co-semigroups together with Riesz 
projections, to analyze the action of the semigroup {e a G( e >} a>0 on vectors of the atomic subspace. This 



analysis leads to a version of Corollary |4.5| for the full microscopic dynamics and Theorem 3.3 follows. □ 



4. Effective atomic dynamics 
In this Section we study the behavior of the solutions of the effective atomic master equation, that is, the initial 



value problem (3.15). Since an important issue of laser technology is to obtain an optical pumping of atomic 
energy levels, we want to understand the time behavior of its solution in relation with the phenomcnological 
Pauli master equation found in standard textbooks on lasers. The small (order A 2 zi7 _1 ) fast oscillations of the 
populations due to the cosine in the master equation preclude a perfect steady behavior at large times. It is 
thus convenient to remove them by averaging the non-autonomous dynamics over a moving period. This leads 
to a pre-master equation proven in Theorem |4.8| whereas the Pauli master equation only extracts the limiting 
behavior of populations at large times (see Section [5]). 

4.1 Methodology and Numerical Illustrations 

To this end, we first represent this non-autonomous evolution as an autonomous dynamics on an enlarged, 
infinite dimensional Hilbert space 

Scvo 2 Sat = B(C d ) . 



See (4.5| below. The latter emerges through an additional degree of freedom which is a new time variable 
denoted by a > 0. 

By iterating the "variation of constants formula", i.e., the integral equation 

Vs,t G R, t>s: f&" = ffj^ + f t^M^t^&v (4.1) 



with the 2irzu 1 -periodic operator 

W (X,V) £ (A lJ7 ) _ £ (0,0) = 7]Cos{mt)&p + A 2 £k g q ( ^ at) _ (4 2) 

we get a representation of f J* as an absolutely convergent (Dyson) series which shows that (t, s) H> is 
continuous and 

Vfc G Z, s, i G K, t > s : T<?f = r^- W ^k ■ ( 4 - 3 ) 

In other words, the dynamics betw een t imes (a + 27rn7 _1 fc) and (t + a + 2i:vj~ 1 k) does not depend on k £ Z. As 
T^lg does not affect populations, (4.2) also shows that they do not change much within a period of the pump. 



As already explained, it is thus natural to average the non-autonomous dynamics over a moving period of length 
2irzu~ 1 to extract the leading dynamical behavior of populations, in particular the inversion of population. The 
latter is described in Section l4~2l 



This first step of the analysis of the solutions of the atomic master equation is quite useful because it enables 
the analysis of the non-autonomous dynamics via an associated (evolution) semigroup denoted by {T a } a >o 
corresponding to an autonomous dynamics. Observe that {T a } a >o acts on an infinite dimensional Hilbert space 
Sevo 2 Sat although the initial non-autonomous dynamics was finite dimensional. Nevertheless, since the 
Hilbert space 23 C .ftat is an invariant subspace of {7~ a }a>o when A = 77 = 0, Kato's perturbation theory |26j 
shows, for sufficiently small coupling constants |A| and |?7|, the existence of a finite dimensional invariant Hilbert 
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space f)Q X ' v ^ of {T a } a >o almost parallel to D. In particular, concerning the dynamics of populations, we can 
finally pass to an autonomous finite dimensional dynamics. Note however that dimS)^'^ > dim 2), as 2) is a 
subspace of the larger invariant space of {T a } a >o when A = r) = 0. This second step is performed in 

Sections lOrEUl 



As a first application, these results are then used at the end of Section |4.4| to study the large time behavior 
of (p(t)). We show in particular the existence of a density matrix p x (Theorem 4.6) which is uniquely 
determined by a balance condition (Theorem 4.7| and approximates P3 (p(t)) properly in the limit t — > 00 
for small enough coupling constants |A|, \rj\. In Section 4.5 we derive an integro-differential equation (pre- 
master equation) on the subspace of block-diagonal density matrices 2) C -f) a t effectively describing the physical 
evolution of the populations. The dynamics of populations is properly described by an integro-differential 
equation and not by an effective differential equation like a Pauli equation. This reflects the strict inequality 
dim^ A,n) >dimS). 

The rigorous proofs can sometimes be tedious and require a number of definitions and notations. Therefore, 
we outline our study by giving now some numerical illustrations of these three master equations, that is: 



(i) The effective atomic master equation (3.151 



(ii) The pre-master equation, see Theorem 4.8 



(iii) The phenomenological Pauli master equation (5.5) as explained in standard textbooks on lasers. See 
Section [5] for more details. 



As explained in Sections [T] [2j an important step of laser technology is to obtain an optical pumping of atomic- 
energy levels. So, we focus here on the so-called inversion of population^] As described in textbooks on 
laser physics, optical pumping is in many situations based on three- or four-level atoms |13j . the second case 
being the more efficient of both. Therefore, we restrict our study on a non-degenerated four-level atom, i.e., 
d = N = 4. In this case, the effective atomic master equation (3.15) is a non-autonomous evolution equation 
on a 16-dimensional Hilbert space which can easily be treated by standard numerical methods. One can then 
understand better the different approximations performed in this section which lead to the pre-master and 
Pauli master equations. 

In our example, the atomic Hamiltonian depends on the parameter w and equals 
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VD J 



(4.4) 



The atom-reservoir interaction is fixed by m = 1, the self-adjoint matrix 



and the coupling function f ± defined on ]R 3 by 

A (p) = H\p\) 
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2tt|p| 



exp(-|p|/2) 



See Section |2.5| The optical pump is modeled here by the time-periodic perturbation 

rjcos(wa) 



/ 1 \ 





V 1 ) 



t e 



3 Note that analytical studies concerning the inversion of population can easily be done, at least for d 
balance equation | |4.45[ l. 



N = 4, by using the 
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of the atomic Hamiltonian H at , see Section 2.3 The initial state of the atom p at eDc f) a t is the Gibbs state 
p g = Pjd (p B ) (2.5) taken at the same inverse temperature f3 as the fermionic reservoir. We set the frequency 
vj, the inverse temperature (3 and the coupling constants A, rj respectively equal to vj = 3, (3 = 0.5, A ~ 0.385, 



and r\ = A 



0.148. 



by 



The population density of the fcth atomic level, at any fixed time a € and for any k G {1, . . . , 4}, is given 
d fc (a) := Pk (p (a)) = Tr C 4(l [H at = E k ] p (a) 1 [H at = E k }) > . 



Here, p(a) is the solution of the effective atomic master equation (3.15) and dfc(a) is directly related to the 
atomic populations because of Theorem 3.3 Similarly, populations can be defined for the solutions of the 
pre-master equation and the Pauli master equation. The plots of all these densities are given in figure [T] and 
clearly show a stable inversion of population as 

d 3 (a) > di(a) > d 4 (a) > d 2 (a) 

for large enough a < 180. As soon as the optical pump is turned off, the systems relaxes to the Gibbs state p fl . 




Figure 1: Illustration of the populations as functions of a = t e at (3 = 0.5, A ~ 0.385, r\ = A 2 ~ 0.148 
and vj — 3. Blue, green, orange, and red lines correspond to the populations of the 1st, the 2nd, the 3rd and 
the 4th atomic energy levels, respectively. The dashed dotted line marks the time a — 180 when the pump is 
turned off, i.e., r\ — for all times beyond this line. The four dashed lines mark the populations of the thermal 
equilibrium state, i.e., the Gibbs state p B at inverse temperature (3 = 0.5 of the thermal reservoir. The two 
figures on the top are computed from the master equation. The pre-master equation gives exactly the same 
picture at this time scale. The figure on the bottom is computed from the Pauli master equation. 

As we can observe in figures [T] and [2j the qualitative difference in the behavior of the solutions of the master, 
pre-master and Pauli master equations for population densities are clear: The pre-master equation removes the 
small, but fast oscillations due to the cosine in the master equation (figure [2| , and the Pauli master equation 
additionally cancels the Rabi (moderate) oscillations still present in the pre-master equation. (Compare the 
plots of figure [T]) . 

Note that Rabi oscillations, which depends on the coupling constant 77, are progressively suppressed by 
dissipative processes. Removing the atom-reservoir interaction by setting A = while keeping 77 ~ 0.148 we 
observe non-suppressed oscillations, see figure [3] In this case the stable inversion of population also disappears, 
as explained in the beginning of Section [2j 

To conclude, our results permit purely quantum mechanical detailed studies of many-level optically active 
impurities used to produce lasing materials. For instance, they can be used to analyze the influence of tem- 
perature and other relevant physical parameters on inversion of population: At least in the example considered 
here, the inversion of population is decreasing with the temperature, i.e., is an increasing function of f3 as shown 
figure |4j and the most efficient parameter to strengthen it is curiously the frequency w and not the parameter 77, 
provided the last is neither too small nor too large. The effect of the degeneracy and the dynamics for very low 
frequencies can also be analyzed. The latter corresponds to a situation close to the adiabatic limit and shows 
an unusual behavior. We postpone this kind of studies to a further paper and start now the rigorous proofs 
related to these numerical observations. 
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Figure 2: Illustration of the populations as functions of a = t £ [0, 2] at f3 — 0.5, A ~ 0.385, r\ — A 2 ~ 0.148 and 
zu = 3. Blue, green, orange, and red lines correspond to the populations of the 1st, the 2nd, the 3rd and the 
4th atomic energy levels, respectively. They are computed from the master equation. The dashed dotted lines 
are the same objects computed from the pre-master equation. 




30(1 



Figure 3: Illustration of the populations as functions of a = t € [0, 300] at /3 = 0.5, A = 0, rj ~ 0.148 and hj = 3. 
Blue, green, orange, and red lines correspond to the populations of the 1st, the 2nd, the 3rd and the 4th atomic 
energy levels, respectively. The dashed dotted line marks the time a = 180 when the pump is turned off, i.e., 
i] = for all times beyond this line. The four dashed line mark the populations of the Gibbs state p g at inverse 
temperature /3 = 0.5 of the thermal reservoir. The figure is computed from the master equation. Note that, for 
a > 180, the system do not relax to the Gibbs state p g at inverse temperature j3 = 0.5 of the thermal reservoir 
because there is no atom-reservoir interaction as A = 0. 
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Figure 4: Illustration of the inversion of population d3(oo) — d.2(oo) for large times, computed from the pre- 
master equation, as a function of the inverse temperature f3 £ [0.2, 20] at A ~ 0.385. In the left and right plots, 
the red lines correspond to r\ = A 2 ~ 0.148, w = 3, whereas the green line represents the choices r\ = A 2 , w — 6 
(left plot) and r\ = 10A 2 , w = 3 (right plot). Note that the atomic Hamiltonian H at depends on the frequency 
w, see (14.41) . 



4.2 From the non autonomous master equation to an autonomous effective dynamics, evolution 
semigroups 

Consider the Hilbert space 

f) cvo := L 2 ([O^irm- 1 ) ,<o at ) (4.5) 
of time-dependent 27rw~ 1 -periodic f) a t-valued functions. The scalar product on $) cvo is naturally defined by 

777 / z& 777 / -ccf 

(/. .9>cvo := ^ J </ (*) . 5 (*)>at dt = to J Trcd ( / (<) * 3 (<) ) dt 



for all f,g £ Sj evo , see (2.6). ,f) a t £ i5evo is seen as the subspace of constant functions of .fjcvo- 

From the continuous two-parameter family {'H^,' , ^}t> s corresponding to the non-autonomous master equation 
(3.15) we can uniquely define a strongly continuous one-parameter semigroup {Ta} a >o by the condition 

Vt G [0, 27TTI7- 1 ) a.e. : T a (f) (t) = f^J (t - a) (4.6) 



for all a > and / £ Sj C vo- Because of ( |4.3[ ), 7^ is an operator acting on f) evo for any a > 0. The strong 
continuity of a 7^ follows from the continuity of (s, t) > r^f 1 ^, and the semigroup property of {7^} Q >o from 
the cocycle property of the two-parameter family {t^^^^s- Moreover, by the norm boundedness (3.16) of 
the evolution family {t[ A s '' 7 ' ) }(> s , the semigroup {T a } a >a C B(Sj cvo ) is uniformly norm bounded: 

VA,?7 G K, a > : || % \\< C , (4.7) 

for some finite constant C £ (0,oo) not depending on A, 77, and a. 

The generator of the strongly continuous semigroup {T a } a >o is the closed unbounded operator 

+ , (4.8) 

sometimes referred to as Rowland operator of the non-autonomous atomic dynamics, with dense domain 

{00 00 \ 

J2 a k e lk ^ : a k £ £ at , £ ||fca fc ||* t < 00 \ C io cvo . (4.9) 
k— — oo k— — oo ) 

Here, 
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in the L 2 ([0, 2ttw : ) , ijat) sense for all 



/ = a k e lk ^ G 2?(G (A '^) , 



and &ivo^ G B(Sj cvo ) is the bounded operator defined, for all / G -f) ovo , by 

Vt g [O^" 1 ) a.e. : £&>>(/) (t) := -C^ (/(*)) 



(4.10) 



Remark 4.1 (The uncoupled atom reservoir case) 

The spectrum ofG^ ' ^ is purely discrete. Therefore, Kato's perturbation theory [26] of discrete eigenvalues can 
be used to study the spectral properties of the generator G 1 --*^ for small \X\ and \rj\. 



The time-behavior of the solution p(t) G Sj at of the non-autonomous master equation (3.15) can be studied 
on the subspace D c B(C d ) = $j at of block-diagonal matrices (cf. (3.1 1) by using the Co-semigroup {T a } a >o- 
More precisely, we prove in the following lemma that, for any (block-diagonal) A G 55, the scalar products of 
the form 

(T a (p),A) cm 

properly describe (i.e. up to small oscillations uniformly bounded in time) the time evolution of 

(p(a),A) at :=Tx(p(a)A) 
whenever the pump frequency vj is sufficiently large or the couplings |A|, |7j| are sufficiently small. 

Lemma 4.1 (Average dynamics over a moving period of length 2ttg7~ 1 ) 

Assume that p at G 3D. The unique solution {p(t)}t>o of the effective atomic master equation (3.15) satisfies the 
bound 



\(p(a),A) at - <T« (p at M)evol < CX 2 ^ 1 1| A|| 



(4.11) 



for all (block-diagonal matrices) A G 3D, all |A| < 1 and all a G Kg". Here, C G (0, oo) is a finite constant not 
depending on p at , A, X, n, vj , and a. 



Proof. By (4.6), 



(p(a),A) at -(T a (p at ),A) t 



w 

2^T 



((€f-r ( t ^l)(pJ,A) dt 

\ /at 



(4.12) 



Therefore, we need to estimate the integrand in this last equality for any a G M.^ and t 6 [0, 2nzu 1 ). To this 
end, we choose, for any a G Rjj" , r (a) G 27ra7~ 1 No such that 

< r (a) - a < 2-kw^ 1 . 
Using the 27nzJ _1 -periodicity (4.3) of the evolution family {r^f^t^s we obtain that 



a,0 



t,t — a ~ 1 ct.O 



r <5+a,<5 



(4.13) 



for any t G [0, 2ttw 1 ] with 

8 := t + r (a) - a G [0, 47rra7~ 1 ] . 
Using the cocycle property of the two-parameter family {T[ X ^}t> s together with (4.13) we have 



JX, n ) /-.(A,»7) 



.(A,r,) w(A,r,) 



t,t-a 



\-S,S 



TO" - 1) + (l-*£3 J 



(4.14) 



with 8 G [0,47ra7 _1 ]. Note that ||[i? at , 
the norms 



l T <5,0 



0{w) [vj = En — Ei G cr([H a t, ■ ])) and thus, we cannot expect 
-1|| and ||l-r£#J| 
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to be small for large pump frequencies w » 1 when 5 — 0(m 1 ). 



However, the integral equation (4.1) implies 



T s+5,s 1 — T S L 



s+S 



(0,0) 9n (A, )7 )A(A,n) tl ,. 



for all set and <5 e [0,47re7 _1 ]. Meanwhile, for all Ae D, 

(i-4°' o) r(A)=o 



(4.15) 



(4.16) 



as (f^ 0,0 -*)* = e 5£at and £ a t(25) = {0}- Recall that the pump is moderate with respect to the atom-reservoir 
interaction, i.e., \rj\ < CX 2 for some fixed constant C € (0,oo), see Assumption 111 He nce, since the evolution 
family {^'^los is uniformly norm bounded (cf. (3.16)), we deduce from (4.2), (4.15) and (4.16) that 



{1-^L)*{A)\\<CX^\\A\ 



for some constant C G (0,oo) not depending on A G 33, A, 77, w, a and S € [0,Attu7 1 ]. By (3.16), 



< C\ 2 m- X \\A\ 



(4.17) 



(4.18) 



for some constant C € (0, 00) not depending on p at , A, A € [—1, 1], 77, 137, a and S G [0, 47rn7 1 ]. 
Similarly, as p at € 25 and ||p at || = 1, 



IK 



r <5,0 



l)(p a t)ll < CA 2 ^- 1 



for some constant C £ (0, 00) not depending on p at , A e [—1,1], 77, ro, and <5 € [0,47rn7 1 ]. Using this together 

(4.19) 



with (3.16) we conclude that 



/ at 

for some constant C € (0, 00) not depending on p at G 25, A G 23, A e [—1, 1], 77, ro, a and 5 € [0, 47rn7 -1 ]. 



From (4.14), (4.18) and (4.19) we obtain 



Vt6 [0.27TW- 1 ] : ((riy-f^l)(p at ),A) < CA 2 ^ 1 ||A|| 

\ /at 

with some constant C E (0, 00) not depending on p at e 25, A G 23, A e [—1, 1], 77, zu, t and a. Combining this 



with (4.12), estimate (4.11) follows 



□ 



Remark 4.2 (General atomic initial states) 

If the density matrix p at € f} at of the initial state uj at is not block-diagonal, i.e., p at G fj at \23, then the assertion 
of Lemma J^.l holds at large times. Indeed, the transient behavior of {p at (i)} tgK + strongly depends on the 



quantum correlations of the initial atomic state, whereas its long time behavior does not depend on the initial 
conditions. The following bound can be shown for arbitrary density matrices p at G -fj a t (i.e., p at G 25 is not 
assumed): 

\(p(a),A) at - (To (pJ,A) e J < C{X 2 w~ l + e— ) \\A\\ 

for some constants c,C G (0,oo) not depending on p at , A, A, rj, w , and a. We omit the details and illustrate 
this fact with a numerical example: see figure^ Indeed, we focus on the study of the dynamics of populations 
and the initial state is in most cases of interest the Gibbs state Q at G 25. For sake of technical simplicity we 
assume that p at G 25 and only observe at this point that all results below on the dynamics {p at (i)} t6R + stay 
correct for all p at G £)at\35 up to a transient factor decaying as e~ ca for some c > 0. 
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Figure 5: Illustration of the populations computed from the master equation as functions of a = t G Kq at 
j3 = 0.5, A ~ 0.385, r\ = A 2 ~ 0.148 and w = 3. Blue, green, orange, and red lines correspond to the populations 
of the 1st, the 2nd, the 3rd and the 4th atomic energy levels, respectively. The dashed dotted line marks the 
time a = 100 when the pump is turned off, i.e., r\ = for all times beyond this line. The four dashed lines are 
the populations computed from the pre-master equation. The figure is computed with initial matrix coefficients 
(Patkfe = 0.5 for (j,k) G {(1, 1), (4,4), (1,4), (4, 1)} and (p at ) jtk = otherwise. I.e., Pat $ D. 



4.3 Dimensional restriction of {T a } a >o 

As explained at the beginning of Section |4j the semigroup {7~ a } a >o acts on an infinite dimensional Hilbert 
space -f)cvo, but the initial conditions we are interested in are constant functions, i.e., elements of the finite 
dimensional subspace i} a t C -f) C vo- It turns out that .f) a t is contained in some finite dimensional subspace -ft:™ 



The latter is invariant with respect to 



which is almost parallel to a finite dimensional subspace Sj-^^ C Sj e 
{T a } a >Q, see (4.26) below. As this semigroup is bounded (cf. ( 4.7 ) ) , the restriction of the autonomous dynamics 
to this invariant subspace describes - up to small errors - the evolution of the solution {p(t)} t >o of the effective 
atomic master equation (3.15). 

To define Sji^v 1 ^ precisely we need some preliminary definitions. We denote by 



1 

2~7ri 



(4.20) 



l*+«l=f 



the Riesz projection [2"6l, Chapter II] associated with the generator G^ X,TI ^ defined by (4.8). Here, 

J?:=min{|e-e'| : e, e' G a([H at , • ]), e ^ e'} > 



(4.21) 



and we assume that the atom-reservoir coupling |A| - and thus |7y|, by Assumption [T] - is sufficiently small to 
ensure that the Riesz projection p^ X '^ is well-defined and has the same dimension as pjP ,0 \ Then, for each 
e G a([H aU ■ ]), 

f.(A,?7) — p(^,V)c, 



is an invariant, 
{7q}q>o- Consequently, 



finite 



dimensional subspace of 



the 



evolution 



Si 



(A,7j) 



= span 




(4.22) 
semigroup 

(4.23) 



is a finite dimensional invariant subspace of {T a } a >o and $)&t C $) 



| n °; 0) and P e ^) 



P, 



(0,0) 



<D{\ 2 ). 
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As a consequence, if |A| is sufficiently small and vo is large enough then the restriction of the semigroup 
{T a }a>o to its invariant space -S^y accurately describes the time evolution of {(p(a), A) at } Q > : 

Lemma 4.2 (Finite dimensional effective autonomous dynamics) 



Assume that p at 6 D. The unique solution {p(t)}t>a of the effective atomic master equation (3.15) satisfies the 
bound 



<P(«M>at " E ( ex P (aG^">Fj A '">) P £ (A '"Vat^ 

ee<T([J? at ,-]) 



/or any (block- diagonal) A G 5), |A| << 1 and a € Rj}~. Here, C € (0, oo) is a finite constant not depending on 
Pat, A > ^, V, ™, and a. 



Proof. Consider the projection 



P, 



onto the invariant subspace Sj^^ ■ Note that fj at C io-n^ an d 



(V)) . = p(A,r;) 

inv ■ / j e 

ecr([H at) .]) 
(0,0) 



d(0,0) 

^inv P = P 



(4.24) 



(4.25) 



By Assumption [l] and Kato's perturbation theory for discrete eigenvalues, there is a constant C € (0, oo) 
not depending on coupling constants A and ry such that 



p(A,n) _ p (o,o) 

inv inv 



(4.26) 



at small |A|. By (4.7), we also observe that the operator family 

is a bounded semigroup for any eigenvalue e_€ er( [H at , • ]) as it is a restriction of the bounded semigroup 
{Ta} a >o onto the invariant subspace Sj{ X ' v ^ (4.22 1. Indeed, 



Va > 0, e € tr([fr at , ■ ]) : exp (aG^P^) = TaP^'" 3 + (1 - P^ ] ). 



(4.27) 



Hence, by (4.24), (4.25) and (4.27), for any p £ £ at , 



T a (p) =T a [{l- P&'*>J Pat ) + T a [P^'" Pl 

= T a (( PT P& n) ) p)+ E ex P (« G^P^) P^p 



>(-M). 



te<7([i?at,-]) 



which, combined with (4.7) and (4.26), in turn implies 



(Ta (Pat) , A) cm - E ( eX P (« G^P^) Pi A '") Pat , A) 
e ee r([iWD 



<G||A||A 2 



with G € (0, oo) not depending on p at , A, A, ?y, n7, and a. The assertion follows now from Lemma 4.1 □ 



As a consequence, we can restrict the autonomous dynamics described by the evolution semigroup {Ta} a >a to 
the finite dimensional subspace iojn' 1 '' 1 C Sj e vo- 
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4.4 Effective block diagonal dynamics 



Since Theorem 3.3 only compares the orthogonal projections P© (p at (t)) (3.2 1 and P® (p(t)) of the atomic 
density matrix p at (t) and pit), respectively, we are only interested in the effective block-diagonal dynamics 
defined by {Pa (p (i))}oo- As shown in the following lemma, this quantity is related to the finite dimensional, 
invariant subspace F)q X defined by (4.221 for e = G a([H at , • ]), see also (3.8). 



Lemma 4.3 (Effective block diagonal dynamics I) 



Assume 



that 



D. 



The effective block-diagonal density 



matrix 



{Pa (p{t))}t>o satisfies the bound 

<P S (p(a)), A) at - (exp (aP^GMpM) ^'Vt,^ 

<CX 2 (1 + vj- 1 )\\A\\ 

for any A £ S) at , \X\ sufficiently small and a £ Kj. Here, C £ (0,oo) is a finite constant not depending on p at , 
A, X, r\, vo , and a. 

Proof. The orthogonal projection Pg acts in Jo at and naturally induces an orthogonal projection, again denoted 
by Pd, in the Hilbert space S) evo by: 

V/ G Toevo, t G [0, 2-kw~ x ) : P s (/) (t) := P v (f(t)) . 

In particular, for any A £ fj at , e £ <r([P at , • ]) and a £ Rq", 

(exp (aG^)p(^ p(^) pauP ^ (j4) ^ 

= (jfc (pj^) exp (aG^P^) Pat ) ,A) cm 

and similarly, 

(p(a),Pv(A)) at = (Pv(p(a)),A) at . 
With the last two equalities we use now Lemma |4.2| to obtain the bound 



(Po (p(a)),A) at - (P S (pf^exp (a G( A <">Pj A <">) Pat ) , A 



<CX 2 (1 + m- 1 ) \\A\\ , (4.28) 

where C £ (0,oo) is some constant not depending on p at , A, A, rj, w, and a. For any e G er([P a t, • ]), note that 

P V P^ = P S P e (°'°) + P s (Pi A '"> - p(°'°>) = 5 e , PD + Pa (Pe (A,7?) " Pi°' 0) ) (4-29) 

with <5 £j£ / being the Kronecker delta. Similar to ( |4.26 1, there is a constant C £ (0, oo) not depending on the 
coupling constants A and n such that 



max 
ee<r([iWD 



p(\v) _ p(o,o) 



< CX Z 



(4.30) 



for |A| sufficiently small. From gTF} and pT2^ -( |L30| ) 

(P S (p(a)),A) at - (P„ (exp (aG^Po^^) p at ) 

< CA 2 (1 + w- 1 ) || A|| (4.31) 
with C G (0, cxd) not depending on p at , A, A, 77, n7, and a. Finally, for all A £ f} at , observe that 

[pJcxpLg^P^p^Xa] 



exp 



exp ya 

>(o,0) , 



A,t/) p(-M) 
J n 



/ evo 



P^'exp aG^Pf ) Pat ,A 
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Using (4.7), (4.271, (4.30) and (4.31 1 the assertion follows. 



□ 



The invariant spaces S)q X '^ (4.22) associated with the projectors Pq A ' i? ' ) are, however, not explicit enough 
for practical purposes. Therefore, the next step is to represent the effective block-diagonal dynamics onto the 
explicitly known eigenspace Sj^°\ To this end, we denote the restriction of G^^ onto the space by 



A_(A,rO ._ p(0,0) G (A,r,)p(0,0) 



(4.32) 



Observe that the semigroups generated by A^'^ and p^'^ Q( x >'i) p^-^ are very close to each other for small 
coupling constants |A| and 

Theorem 4.4 (Uniform norm estimates on semigroups) 

For any vo G R, e G (0, 1), and any a G Rq , there is a constant C ro . e G (0,oo) not depending on X, r\, and a 
such that 

exp (aA( A >">) - exp (aP^ G^ P^) || < C^ e |A| 2(1 " £) . 



The proof of Theorem 4.4 needs some technical preparations. For the sake of clarity we defer it to Section 6.2 
We stress that C^e can vary considerably with the pump frequency w. 



Thus, Lemma 4.3 and Theorem 4.4 imply that the evolution of Pjj (p (t)) can be approximated by the 
semigroup generated by A^'^ : 



Corollary 4.5 (Effective block diagonal dynamics II) 

Assume that p at G 2). The effective block-diagonal density matrix 

{Pd (p{t))}t>o satisfies the bound 



(P s - (P ( °< 0) exp («A(^)) P ( °'Vt,A) c 



< CL« I A 



2(1-6) 



\A\\, 



for any A G Sj ai , |A| sufficiently small, any e G (0,1) and a G Rq . Here, C ro E G (0, oo) is a finite constant 
depending on the pump frequency w and e but not on p at7 A, \, rj, and a. 



In Corollary |4.5| we approximate th e effe ctive block-diagonal dynamics by some time-evolution on the finite 
dimensional subspace ^q 0,0 ^ £ fj cvo (4.22). This Hilbert space is not a subspace of f) at £ f) vo because of 
oscillating terms present in it, but it can be explicitly defined as follows. 

Recall that the eigenspaces of the atomic Hamiltonian P at G £>(C d ), associated with the eigenvalues for 
fc G {1, . . . , N}, and their dimensions are denoted by Hk C C d and n k G N, respectively. See Section 2.2 By 
taking any arbitrary orthonormal basis {e^ °f ^fc f° r eacn k € {1,..., N} we define the elements 



(fe'V) 



for any k, k' G {1, . . . , N}, n G {1, . . . , n^} and n' G {1, . . . , ?ifc'} by 
Vfc"G{l,...,7V}, n"e{l,...,n fc »}: 



(kV) (fc") 
(fc,n) n" 



(*') 



W» ; := <5n,r l "^fc,fe"e„/ 



(4.33) 



Then, straightforward computations show that the Hilbert space ^q ' ^ equals 



fjo°'°^ = span 



j e -i*(*V- 



>w, 



(k',n') 
(k,n) 



(fc, fc') G t_ ro U to U in, n G {1, . . . , n k }, 



^' G {l,...,^}}. 



Recall that t mro is the set defined by (3.7) for e = mw G a([H &t , ■ ]) and mG {— 1,0,1}, that is explicitly, 

i-„={(l,N)}, t = {0',i):ie{l,2 > ...,JV}}, t ro = {(7V,l)}, (4.34) 



2G 



because w :— En — E\ > 0, see (2.10). 
Obviously, 

,(o,o) 



4°' 0) nio at = P s (i3at) 



and SJq ^ -fjat- Nevertheless, we can remove the oscillating terms by defining a unitary map Uo from jo 
to the atomic subspace 



(0,0) 



P(0,0) fTJ/( fc '^ 1 ') 

f)n .= span<| WA. „^ y 



\ (*>■ 



j y |(/c,fc')et- CT Ut Ut ro , ne {l,...,n k }, 



/€{l,...,n fe -}} C^at (4.35) 



as follows: 



U n e" 



(fc,n) 



(4.36) 



for any k, k' G {1, . . . , N}, n G {1, . . . , n k } and n' G {1, . . . , n fc .}- Clearly, P s (ij at ) C ^ °' 0) as U Pd (.flat) 
Pd (Sat)- 



Hence, by Corollary 4.5 the behavior of Po (p (a)) can be studied through the semigroup acting on fig 0,0 " 1 C 
flat and generated by the operator 

A( A >") := UoA^^Uq- . (4.37) 
This follows, for any initial density matrix p at G .flat, from the equality 

P s (exp (aA^) P (0 ' 0) (p at )) = P 3 (U^exp (aA( A >">) U P S (p at )) 



P 



(exp (<*A< A -">) P S (p at ) 



(4.38) 



As a first application of the above results we are now in position to study the large time behavior of the 
effective block-diagonal density matrix p£> (p (t)): 



Theorem 4.6 (Large time behavior of Pg {p(t))) 
(i) For all p G fj °' 0) 



Jo ' 



A( A >") (p)^£ p (p) + A 2 £ TC (p)Gf, ( 



(0,0) 



wi£/i £ p and Ztz defined by (3.5) and (3.6), respectively. 

(ii) Assume that p at G T>. There is a unique density matrix G fj '°^ such that JS^'V) (p^) = and 



limsup ||P S (p(a)) - Po foJH < |A| 



2(l-e) 



(4.39) 



/or all |A| sufficiently small and any e G (0, 1). Here, C ro . £ G (0, oo) is a finite constant depending on zu and e 
but not on p at , X, rj. 



Proof, (i) For any p G i^o' ^ an d t G [0, 2ir-cu *) a.e., note that 

d 



dt 



[XJ* oP ] (i) + £ at \U* oP ] (t) = 



Then, by construction, 



G &«)pMv*p] (t) = \ (c rat + e- rat ) £ p [U2p] (i) + A 2 £ TC [U$p] (t) 



(4.40) 



for any p G fj ' J and i G [O^Trcc" 1 ) a.e., see (3.3), (4.8 )— (4.10 ) and (4.20 1-( 4.22 1. Using the explicit expressions 
(3.12) and (3.13) one easily checks [£r,,Uo] = 0, which obviously implies that 

Uo£r,Uq = £7?, • 
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On the other hand, we have 



(1 - 1 [£ p = 0]) [XJ* p] (t) = 



P\M e 


3- u ™Pn.i (p) P 




where 1 [£ p = 0] is the projection onto the kernel of £ p and Pj t k are the orthogonal projections 



P 



J,* 



Ej] 1 [flat 



for all j, k € {1, ... , ./V}. Hence, we obtain 



JSfc] G B(Sat) 



(4.41) 



P 



(0,0) / jtro 



(I 



£ P [Uop] (*) 



Pi.jv (p) hp - h*P Ntl (p) 

Atw (Pn,n (p) ftp - W (p)) 



e ltro (P,i (p) ft; - ft;Pv,iv (P)) N 


Pjv,i (p) h* ~ ftpPi.jv (p) 



Multiplying now Uo from the left in (4.40) and using the computations above we arrive at 

Uo p(o.o) G r(A,n)p(o ) o) USp = ^ (p) + ^ (p) 



for any p G ^q 0,0 '. 

(ii) Note that the bounded operator A( A,I? ' given by the equality of the first assertion (i) makes sense for all 
ato mic density matrices p € jj at , whereas i^Q 0,0 ' C fj a t is an invariant space of A^-''' g i3(ij at ). By Remark 



3.1 



Assumption [5J Theorem 



6.10 



and assertion (i), A^'^ 6 B(Sj a t) is the generator of a relaxing Markov CP 
semigroup for all A, r\ <E M 2 , A ^ 0. See Definition 6.9 and Remark 3.2 I.e., there is a unique density matrix 
Poo <= #at such that, for any p € f) at , 



lim 

a— ^oo 



(exp (aA^">) p) = Po 



(4.42) 



6.10 



As Sj 



(0,0) 



It follows that p^ € £j a t is the unique density matrix satisfying A"''" (p^) = 0, see Theorem 
an invariant space of A^'^ € B($3a,t) containing density matrices, one must have p^ € i^ ' ''. Using (4.38) and 
the fact that Po (p at ) G i} a t is also a density matrix, we obtain 



lim P, 



(0,0) 







exp 



(aA<*">) P (0 < 0) (p at ) = P s (p c 



The inequality (4.39) then results from Corollary 4.5 and the finite dimensionality of the Hilbert space Sj at . □ 

It now remains to characterize more precisely the block diagonal projection p x := Pj) (p^) € of the density 
matrix ~p x of Theorem 4.6 (ii). This is done via a balance condition for populations in Theorem 4.7 below. 

We first define the orthogonal projection 



P T 



Ps(o,o) 

'' J 



Pj> € B(9) a 



and the operator 



i.e., £ = P^AO'^P^-. 



<£ := P^ (|£ p + A 2 £ K ) P^ e B(i5 at ) , 



The operator £ = <L\ only depends on A € R: Observe that Pjj and P^ are orthogonal projections onto the 
Hilbert space 5$'^ C i} a t for which ran(Pjj) and ran(P^) are invariant spaces of £7?,. The operator £ p however 
maps the subspace ran(Pg) of block diagonal matrices to the subspace ran(P^) of off-diagonal matrices ini^ 0,0 -' 
and vice versa. In particular, 

£ p Po = P^£ P P S and P£,£ p Ps, = 0, (4.43) 
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and hence, 



€ = € x = \ 2 P££, R P£ and <t 



A 2 P^ £^P^ 



for all A,t)Gl. 

Recall that, by Remarks 3.1 3.2| Assumption [2] and Theorem 6.10 the atom-reservoir Lindbladian Z-r is 
the generator of a relaxing Markov CP semigroup on .f) at and all non-zero elements p £ ct(£r,)\ {0} have a 
strictly negative real part Re (p) < 0. On the other hand, explicit computations show that the density matrix 
p g = P<z (p g ) of the Gibbs state Q at belongs to the kernel of £-r, i-e., &n (p s ) = 0, provided that the parameter 
/3 G (0, oo) in (|2.5[) is the inverse temperature of the reservoir. Thus, 



ker £ K = 



C D and ker Z* n = C • l € d C £ . 



Note that the second equality is an obvious consequence of Theorem 6.7 



It follows that the operator £ and its adjoint €* are both invertible on the subspace ran(P^ ). Therefore, 
using (4.43) and standard results on Feshbach maps [37J Theorem 2.1, Remark 2.2], we deduce that p M- P-s{p) 
defines a one-to-one map from 



to the space 



with 



ker A^> n ^°< 0) = ker (a 2 £ k + ^£ p ) D ^ 



ker ( X £,iz + ^ <B p ,tj ) i ? 



(0,0) 



(4.44) 



By uniqueness of the density matrix £ ^q 0,0 " 1 satisfying A^' 17 ^^) = (Theorem 4.6 (ii)), the following 
unique characterization of the populations — Pg (p^) holds: 

Theorem 4.7 (Characterization of p^ via a balance condition) 

Poq = Pj) (poo) «s tte unique (block diagonal) density matrix p^ € 2) satisfying the balance condition 



£k (Poo) + ^4 S P. K (/>° 







(4.45) 



As explained above, all non-zero elements p G ct(£r,)\ {0} have a strictly negative real part Re (p) < 0. In 
particular, there are constants C, c € (0, oo) such that 



PiVH < Ce" sc 



(4.46) 



As a consequence, by using (4.43) and expressing resolvents of generator of semigroups through Laplace trans- 
form, we can rewrite the operator 25 Pi 7?., which describes the pump contribution to the (quasi-) steady popula- 
tions, as 



/•OO 

®p,tc - / £ P e s£K £ p P D ds g B(«at) ■ 
Jo 

This formulation is important to get the balance condition from a dynamical principle. 



(4.47) 



4.5 The pre master equation and the balance condition 

We analyze in the present section the time-evolution of density matrices 



Since 



p s (a) := P S exp (aA^) P s (pj € P s (^°' 0) ) 
\\Pv(p(a))-p 3 (a)\\<C^ £ \X\ 2(1 - e) 



(4.48) 
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with C ro , e € (0,oo), e € (0,1), not depending on p at , A, A, 77, and a (cf. Corollary 4.5 and (4.38)), the 



density matrix (a) accurately describes the real evolution of atomic populations at small couplings. One 



important consequence of this analysis is the derivation, after Theorem 4.8 of the balance condition (4.45 1 from 



a dynamical principle. Indeed, p^ (a) satisfies an integro-differential equation, called pre-master equation, and 



Theorem 4.7 is equivalent to the fact that p^ is the unique stationary state of the Markov approximation of 
this integro-differential equation. 

Theorem 4.8 (The pre master equation) 

The family {p^ (a)} a >o <Z T) of block-diagonal density matrices obeys the integro-differential equation 



2 

Va>0: (a) = A 2 £ K (p s («)) + ^2 / V s£tc £ p (p s (a - sA" 2 )) ds 



(4.49) 



with p s (0) = p at e ID. 

Proof. The proof of this assertion is standard (see for instance [38J Chapter 7]) and is given here for completeness. 
The two-fold iteration of "variation of constants formula" yields the equality 



ds\ 



ds 



(4.50) 



ds 2 e (Q " Sl)A £TC £ n e (si " S2)A £w £ n e S2 



Using that 
and 



wPb&n c p£&r\ )) c p^°< o) ) 



we readily deduce from (4.50) that 



dai / ds 2 ef-^^^^et^^^^^P^^'^P; 



Deriving this last equation we get 
d 



da 



Pv (a) = A 2 £ K fa (a)) + ^ J £ p c sA2£ -£ p (p s (a - s)) As 



from which we deduce the theorem by a trivial change of variable. 



□ 



By combining (4.42) with the equality A^ 1 ') (p^) = (Theorem 4.6 (ii)), the density matrix p^ (a) must 



converge to P35 (p^) and its derivative must vanish in the limit a — >• 00. By (4.47) and Lebesgue's dominated 



convergence theorem, the limit — P55 (fi^) must solve the balance condition (4.45), as already proven in 
Theorem EH 



Remark 4.3 (Moderate optical pump) 

The balance condition shows that the contribution of the pump to the final atomic state is of order rj 2 / \ whereas 
the contribution of the atom-reservoir interaction is of order one (in the parameter rj 2 /\ ). As explained in 
2.1, this justifies Assumption [7J that is, \r/\ < CA 2 . In particular, this regime follows Observation (b) 



Section 



given at the beginning of Section^ 



5. Generalized Einstein Coefficients and Pauli equations 
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Because of Theorem 4.7 we can interpret A 2 21k, with 



2l K := Ps^k-Pd , 

2 

as spontaneous transition rates and jW5B p ,tc as effective stimulated rates between atomic energy levels. In 
order to make this precise it is natural to impose that and 2$ Pj -r, generate Markov semigroups on T) which 
preserves positivity. 

The operator 21k has this property for any choice of parameters because £-r, generates a Markov CP semigroup, 
which preserves the subspace 2). Note however that this feature is in general not satisfied by the operator *B P; tj,. 
A simple counter-example with N = 2 and d = 3, where one energy level is two-fold degenerated, is given in 
Section |6.4| This fact is not very surprising. Indeed, as discussed in Section |4.5[ the balance condition comes 
from a Markov approximation of the restriction of a CP dynamics. It is well-known that this kind of construction 
can destroy positivity [31 Section III.l]. 

Hence, we will assume in this section the following: 

Assumption 3 (5B Pi k as transition rates) 



The operator 25 Pj -r, defined by (4-44) * s the restriction on 2) of the generator of a Markov CP semigroup on 
B(C d ) with an invariant space T) dB(C d ). In particular, S Pi k generates a Markov semigroup on D which 
preserves positivity. 



A sufficient condition on the Lindbladian Zn to satisfy Assumption [3] is given by Theorem 6.11 in Section 6.4 



It is always satisfied if the 1st and the Ath atomic energy levels are non-degenerate. Indeed, the condition 
stated in Theorem |6 . 1 1| physically corresponds to the following: 

• The pump is uniformly resonant, i.e., the reservoir-impurity interaction does not split the spectral line 
corresponding to the A-l atomic transition. This atomic spectral line may however move as a whole 
under the influence of the reservoir (uniform Lamb shift). 

• The decoherence time is uniform for the A-l correlations, i.e., the reservoir does not induce a splitting of 
ran(Pjv,i) in smaller independent coherence subspaces. 



Assuming from now Assumption[3j we are in position to define in Section 5.1 what we call generalized Einstein 



coefficients. These coefficients yield the Pauli master equation and Einstein's relations, respectively described 
in Sections 15.21 and 15.31 



5.1 Generalized Einstein coefficients 

Recall that the left and right multiplications are 
and the orthogonal projections Pj ^ are defined, for all j, k € {1, . . . , A}, by (4.41 1. Clearly, 



defined 



by (2.3), 



N 
J,fc=l 



and 



2l K = 2^ A J-k 

j,k£{l,...,N} 



(5.1) 



where, for all j, k £ {1, . . . , A}, A^j. := Pj.j£-izPk,k- Since generates a semigroup which always preserves 
positivity, for any j, k € {1, . . . , A} such that j ^ fc, the operator Aj^ := X 2 Aj^ defines a map from B + ('Hk) to 
B + {%j ) and hence can be interpreted as the spontaneous transition rate from the kth to the jth atomic energy 
level. Here, S + (f)) denotes the set of positive operators on the Hilbert space fj. The operator Aj j := A 2 A J :) - 
for any j 6 {1, . . . , A} is then responsible for the trace preservation of the total dynamics generated by the 
operator A 2 21-r,. 

The spontaneous transition rates Aj^ can explicitly be computed from the quantities defining the microscopic 



model. By (3.4), (3.6), (3.11), (3.12| and (3.13), we have 



A 3 - k - 2 (1 - 5 jtk ) 



i=i 



2c W yW 



v k,k 



N 

i=i 



l.k 



y(l)*y(t) , yyj'y 
"Ik ' h ' ' h ' 



' l.k 



l.k "l.k 
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for all j, k € {1, ... , N} and with being the Kronecker delta. With this expression, and the corresponding 
ones (5.2) for Bj- j. below, we see what role is played by the dissipative effects due to the atom-reservoir 



interaction for the behavior of (quasi-) steady populations of optically pumped atomic energy levels. This is 
already mentioned in Observation (a) at the beginning of Section [2] 

Analogously, define 

Bjv,! := ~PN,Nh p (pn,1 (P-Ira 

—Pn,nIi^ \ Pi,n 
Bi,jv := -Pi,ih p (Pi tN (£\ ts 



ran(f ^ 
ran(P^ 



-1 



/\ i ) h^P 1A 
P\..\ ) h p Pi,i , 



ran(P^) J P 1,N J h^P N ^ N 

-Pl,lh^ (pn,1 (£| ran (f£)) P N,l) hpPN,N , 

-1 



&N,N :— P/V,Af^p \^Pl,N (^lran(P^) ) Pl,N ) h p PN,N 



+P N , N h* [ (£ ra 



P 



4 

h n P\ 



Bi,i := I P N 



,i (a 



•an(P^) 



f Pi.ijp [Pi,n (c:; rau( 



A . ! j 



'V.\ ) gPl.l 



(5.2) 



Recall that /i p € £?(C d ) maps "Hi to Hat and its kernel equals ker (h p ) = Hi, whereas h* maps Hn to "Hi with 
ker (hp 1 ) = Tij^, see Section 2.3 Observe meanwhile that ran (Pi,at) , ran (Pn,i) C ran(Pg ) are invariant spaces 
of the Lindbladian Zn- As a consequence, from (2.11), (3.5) and (4.44), 



23p,TC — Bjv,i + Bi.at + 'Bn.n + Bi ; i 



(5.3) 



For any j, k £ {1, N} such that j ^ fc, the operator 



maps B + (Hk) to B + (Hj) and is associated with the stimulated transition rate from the kth to the jth atomic 

2 

energy level. Similar to the spontaneous transition rates, the operators Bj j := ^j^Bj j are such that the full 

2 

dynamics generated by j^^p^n preserves traces. 

The objects A^k, P>j,k seen as maps — > Kj between cones Kk C B + CHk), K,j c B + (Hj) are called here 
generalized Einstein coefficients. They satisfy strong constraints, named here generalized Einstein relations, 
which have consequences for the structure of the (quasi-) steady populations p x through the corresponding 
balance condition satisfied by p^. These relations are discussed in Section 5.3 below, after introducing the 
Pauli master equation. 



5.2 The Pauli master equation 

By Theorem |4.7| we can see the unique (quasi-) steady populations p^ as the stationary state of the phe- 
nomenological (quantum) Pauli master equation 



2 



V*>0: ^Q(t) = \ 2 Vl n (Q(t)) + ^® p , n (Q(t)), g(0) = P s (p at ) , (5.4) 
at 4A 

which generalizes the usual (classical) Pauli equation for populations found in standard textbooks on laser 
physics. 

From Assumption [3] and Theorem 6.10 the unique solution 



g(t) = Pxi ( £?(£)) of the Pauli master equation converges to the unique density matrix p x solution of 
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the balance condition (4.45) as t — > oo. By Theorems 3.3 and 4.6 (ii), one thus extracts from the Pauli master 



equation the correct asymptotic behavior of the atomic dynamics at small couplings and large times: 

limsup 1 1 Pa, ( Pat (i)) - g(t)\\ < C CT , e |A| 2(1 - £) 

t— J-oo 

for some finite constant C roj£ G (0, oo) depending on w and e G (0, 1) but not on p at , A, and r\. 

Observe also that the Pauli master equation gives the true effective dynamics when r\ = 0, i.e., when there 
is no optical pump. In this case, the master, pre-master and Pauli equations are the same. However, for 



rj =/= 0, the dynamics governed by (5.4) is generically quite different from the time evolution of the density 
matrix {p v (a)}ct>o, which corresponds - up to small errors - to the real atomic dynamics {P® (|0 a t(i))}t>o> cf- 



Theorem |3.3| and 4.6 (ii). See, for instance, the numerical examples in Section 



4.1 



The reason for such a discrepancy are the memory effects related to the stimulated processes: {p s (a)} a 



>o 



is governed by a (pre-master) integro-differential equation (4.49), the Pauli master equation being its Markov 
approximation. Thus, the smaller the decoherence rates of the N—l correlations as compared to the spontaneous 
transition rates Aj jk of the atom, the less accurate is the dynamics given by the Pauli master equation compared 
to the microscopic atomic dynamics. (The decoherence rate and Lamb shift of the N—l correlations correspond 
6.11| to I Re{£jv,i}| an d I m {CAf,i}: respectively.) 



in Theorem 



By Assumption |3j the solution of ( |5.4[ ) evolves in the positive cone 

S+ := B+(C d ) n 2) = co / (J B+(n k )^ . 

Here, co (m) stands for the convex hull of the set m. By using the spontaneous and stimulated atomic transition 
rates Aj, k ,B jik 
3 e {!,..., N} : 



B + {%k) — > B (Hj) defined in Section 5.1 the Pauli master equation (5.4) reads, for all 



V* > 



d 
df 



37 <?,"(*) 



N 

£ 



(A ijfc + B hk ) g k (t), gJO) - Qj G B+CHj) 



(5.5) 



with gjt) = P jd (g(t)) and where B hk := if {j, k} £ {1, N}. 



In many situations, for instance in presence of symmetries, Aj k and Bj k define maps from subcones K k C 



B + (H k ) to subcones Kj C B + (Hj) and we can write the evolution equation (5.5) with initial conditions 
Qj(0) = Qj E Kj for j G {1, ■ • • ,N}. The unique solution of this initial value problem satisfies in this case 
{Qj(t)}t>o C Kj for all j G {1, . . . , ^V}. If the initial density matrix g at G f) at is chosen such that Pjj (g at ) G Kj 
for all j G {1, ... , N}, the jth population Qj(t) then converges to Pjj (g^) C Kj. This system of N differential 
equations is then the Pauli master equation of the invariant family {K k } k=1 of subcones K k C B + ('H k ). 

A simple, sufficient and necessary condition on subcones {K k }^ =1 to ensure that Aj_ k and Bj k both map K k 
to Kj is given in the following definition: 

Definition 5.1 (Invariant family of cones) 

A family {K k }^ =1 of subcones K k C B + (T-L k ) is an invariant family whenever 

K := co ^(J K^ 



is invariant under the action of the semigroups {e* 3171 } (> and {e t>8p ' R } 



t>o 



By the Trotter product formula, observe that the subset K defined in this definition is also invariant under the 
action of the semigroup 



exp t(\ z % n 



]lf_ 
4A 



t>0 



Furthermore, the invariance of the family {K k } k=1 yields 

Vj, k e {!,..., N}: Aj tk (K k ),Bj tk (K h ) c A , 
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Conversely, if a family {K k } k=1 of subcones K k C B + (H k ) is such that Aj >k and Bj tk map -ftTjt to Kj then 
{-^fcl^Li is clearly an invariant family. 

One trivial example of an invariant family of cones is given by taking K k = B + (Hk) for all k 6 {1, . . . , N}. 
But the smaller the dimension of such cones is, the more classical is the description of the final state via the 
transition rates {Aj^} 1 ^ k=1 and {Bj, k }^ k=1 . It can even happen that the dimension of all subcones {K k } k=1 can 
be chosen to be one. The latter trivially happens with K k = B + (H k ) when all atomic levels are non-degenerated, 
i.e., the dimension n k of the eigenspace T-L k is one for all k G {1, ... , N}. Note however that the non-degeneracy 
of all atomic energy levels is not a necessary condition for the existence of such one-dimensional cones and 
under certain circumstances the fully classical picture of the (quasi-) steady populations p^ is valid. In such 
a case the results of this paper directly relate the coefficients corresponding to the classical Pauli equation to 
microscopic quantities. 



5.3 Einstein's relations 



By (5.2H5.3I 



B i,k = V 2 FjA^N,i,A hN ), with Aj, k := A"/',.,.i2,\> /',.;, . 

In other words, Bj k is proportional to the intensity i] 2 of the pumping (monochromatic) light and proportional 
to a fixed function Fj t k(Afj,i, ^-i,n) of the spontaneous "off-diagonal" transition rates A^^ and A\ jq. Einstein 
derived similar relations, called here Einstein AB-relations, for an atom interacting with a (broad-band, i.e., 
non-monochromatic) black-body radiation field in his seminal paper [11] . This was performed by using phe- 
nomenological considerations about the expected final state of the atomic populations and the asymptotics of 
the light intensity at large wave-numbers (Maxwell distribution). 

Note that Fj, k strongly depends on the specific setting. The function Fj ik appearing in the present paper 
cannot be compared to the one appearing in Einstein's work. However, the fact that the stimulated, operator- 
value coefficients Bj k only depend on light intensity and spontaneous (operator-valued) coefficients Aj k seems 
to be universal. We stress that this property is rigorously derived here from a microscopic quantum mechanical 
description of the system under consideration and not from phenomenological assumptions. 

Einstein also gives in his work a relation between the stimulated transition rates Bj k and B k f Denoting the 
degeneracy of the fcth atomic level by n k , he obtained the equations 

Vj, k € {1, . . . , N} : n k B jM = n J!,., , 

named here Einstein BB-relations. 

Let p k denote the population in the fcth atomic level and define the stimulated flux from the fcth to jth atomic 
level by fj ik := Bj k p k - Then the Einstein BB-relations for fluxes reads 

Vj, k € {1, . . . , iV} : I',,"!. /,.!.- - pmjfkj = . (5.6) 
In our setting, a variant of BB-relations holds for fluxes, at least for density matrices in the subcone 

K°:=co(\jK°) , 



\fe=i 



where, for all k G {1, . . . , N}, 



Kl := K+ • 1 [H at = E k ] C B+{U k ) 



In this context and for any density matrix p £ fj at , the population in the kth atomic level is naturally defined 
to be 

Pk (p) = Tr cd {P Kk (p)) = Tr cd (l [H at = E k ] pi [H at = E k ]) > 



for k € {1, ... , N}, see (|2J]). Similarly, for all j, k € {1, ... , N}, 

fj,k(p) '■= Tr c d ( B j,k{p)) = Tr C d (B jtk P kjk (p)) 

represents the stimulated flux from the fcth to the jth atomic level with respect to the density matrix p £ Sj at - 
Then one proves the following: 
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Proposition 5.2 (Einstein BB relations for states in K°) 

For any p G K° , 

Vj, k G {1, . . . , N} : Pj(p)nkfj,k(p) ~ Pk{p)njf k ,j{p) = 
Proof. Clearly, for any k G {1, . . . , N} and p e K°, 

n k Pk,k (p) = [-Hat = ^fc] • 

As a consequence, it suffices to prove the equality 

TV (B jtk (l [H at = E k })) = Tr cd (B kJ (l [H at = Ej})) . 



(5.7) 



Recall that Bj k equals Bj k = for all j G {2, . . . , TV — 1} or k G {2, . . . , N — 1}. So, we only need to prove 



< pT7| for j, fc G {1, N} and j ^ fc. By Q , 

2 

#l,iV (1 [#at = ^JV]) = — TT 2 

4A 



(/i;(£| 1 . an(P . ) )- 1 (/ lp ) + (ci^px))- 1 ^)^) 



B N ,l (1 [fl. t = = (/*p( C li»»(^)) ^p) + (^IranCP-))" 1 !^)^ 



Therefore, using the cyclicity of the trace and the identity 



A € ran(P^ ) C i? a 



the assertion (5.7) follows for j, k G {1, AT} and j ^ k. 



□ 



Now, without any optical pump, i.e., for r\ = 0, the final density matrix p^ belongs to the subcone K°. More 
precisely, in this case p^ — p g G K is the atomic Gibbs state (2.5 1 with same inverse temperature j3 as the 
one of the reservoir. As a consequence, wc infer from Proposition 5.2 together with Kato's perturbation theory 



[2^] for non-degenerated eigenvectors that approximated Einstein BB-relations hold for the (quasi-) steady 
populations p^ for weak pumps. 

Corollary 5.3 (Einstein BB relations at weak optical pump) 

For all (?7, A) G W 2 , 

rj 2 

Vj,ke{l,...,N}: \Pj{Poo) n kf]APoo)-Pk{Poo) n jfk,j{Poo)\<C^Y ■ 
Here, C ro G (0, oo) is a constant depending on vj but not on j, k, p &t , A and r\. 



6. Appendix 



We give in Section 6.1 a heuristic derivation of the time-dependent Lindbladian 2^' ,r ^ (3.3 1 driving the 
effective atomic master equation ( |3.15 1 as well as the proof of Theorem 4.4 in Section 6.2 Section|6.3|is a short 



review on completely positive (CP) semigroups by focussing on results which are relevant for our analysis in 



order to facilitate the reading of the paper. Finally, in Section 6.4 we discuss positivity questions related to the 



interpretation as transitions rates of the coefficients of the balance condition ( 4.45 ) characterizing uniquely the 
(quasi-) steady populations. 



6.1 Heuristic derivation of the time dependent Lindbladian 

1. It is natural to use the GNS representation (fj, n, Qq) of the initial state uio := cj at 
because it provides a Hilbert space structure which simplifies the analysis 40J. Note that $j := Sj at 
n := 7r at ® n-iz and Qq := Q a t ® fl-R, where ($) a t, 7r a t, ^at) and (Sin, ir-R, ^iz) are the GNS representations of w a t 
and Lo-n- As w a t is assumed to be faithful, the initial state ujq is also faithful. In particular, n is injective. So, for 
simplicity, it (A) and ir (V) are denoted by A and V. Moreover, the cyclic vector Q Q of the GNS representation 



uj-jz (|2.15| of the system 
Sin, 
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is in this case separating for 971, i.e., Atlo = implies A = 0. Here, the von Neumann algebra 97t := V" is 
the weak closure of the C*-algebra tt (V). The state u>o on V extends uniquely to a normal state on 971. This 
extension is again denoted by wo- 
lf w a t = flat is the Gibbs state (cf. (|2.5|) then loq is a (/3, r)-KMS state, where r := {r t } te u is the strongly 



continuous one-parameter group of automorphisms of V defined by (2.14 1. As a consequence, by [25, Corollary 
5.3.4], r also extends uniquely to a a-weakly continuous ^-automorphism group on 971, which is also denoted 
by t :— {T t }t e ffi. Moreover, this extension has a unique unitary representation {Ut}teR, i.e., 

Vf G M, A G 9n : r t (A) = J7 t A[/ t * , 

such that Ut^lo — £Iq. As t i— > r t is cr-weakly continuous, the map t i-> Z7 t is strongly continuous. Therefore, 
the unitary group {Ut}t^R has an anti-self-adjoint operator iL as generator, i.e., C/t = e ltL . In particular, 
O G Dom(L) and L annihilates O , i.e., L£l = 0. 

Now, in the case where w a t is any faithful state of the atom, observe that the space fj := Sj^t <8> fin 
as well as the map it := 7r at (g) do not depend on w at . Only the cyclic vector f2 := ^at ® ^iz of the 
GNS representation of loq changes. In particular, the von Neumann algebra 971 and the corresponding 
extension of the *-automorphism group r do not depend on the initial state w at of the atom. Therefore, 
we can take any faithful state w at and use the operators {T t }tem. acting on 971, and the operators {Ut}tem 
and L acting on S) defined above for w at = g at . Note however that the equality LQ, = does not hold, in general. 

2. We then use the (Tomita-Takesaki) modular objects A, J of the pair (971, Q,q) to generalize [ID], where the 
notion of C-Liouvillcans has been first introduced to investigate non-equilibrium stationary states, to the case 
of time-dependent C-Liouvilleans. For the detailed proofs of the results stated above we refer to the paper [T]. 
See also [5] for further applications of time-dependent C-Liouvilleans. 

The pump term and the atom-reservoir interaction are taken into account by defining a time-dependent 
C-Liouvillean 



with 



:= iL + iH^® 1^ - tJA^ 2 {H^ <8 1^)A-^ 2 J + iW t - UA 1,2 W t A" 1 / 2 J (6.1) 

m 

W t :=r]COs{mt)H p + \J2Q t »*(/<) • 



Indeed, we can prove that the time-dependent C-Liouvillean C t generates an evolution family {U^'^y^a 



representing in the GNS space fj the strongly continuous two-parameter family {r t s ' }t> s generated by the 



t,s 



symmetric derivation s\ X (2.17). More precisely, we show that the non-autonomous evolution equations 



Vt>s: dtUft") = CtU$ v) , d s U^ v) =-U^ v) C s , / . 1 (6.2) 

on Dom (L) have a unique solution {U^ v) } t > s C B($j). The operator u£ v) possesses a bounded inverse 
(uj: s''''') -1 for any t > s and is 2ttzu~ 1 -periodic in the following sense: 

Vi > a : ufr' v) = U^f , „ . 

Then, one finally proves that U^^Qq = Oo an d 

Vt > a, A G V : (A) = EZ^A^)" 1 . 

3. To exploit these facts and obtain the atomic time-dependent Lindbladian described in Section [3T| we use an 
explicit GNS representation of the initial state under consideration. As explained above, this GNS representation 
includes of course the GNS representation of w at already described in Section |2.2| and so, it remains to give 



3G 



the an explicit GNS representation (Sj-jz , ttr , il-jz ) of the (faithful) quasi-free equilibrium state uj-ji (2.2 1 of the 
fcrmionic reservoir at inverse temperature f3 g (0, oo). 

The usual explicit form of the GNS representation of fermion quasi-free states is the Araki-Wyss represen- 
tation |39j . Here we use however the (unitarily equivalent) representation introduced by Jaksic and Pillet (40j . 
the so-called Jaksic-Pillet glued representation, because it is well-adapted to spectral deformation methods. 

So, consider the Hilbert space 

f) 2 := L 2 {R x S 2 ,C) , (6.3) 

where S 2 C M 3 is the two-dimensional unit sphere centered at the origin and K x S 2 (spherical coordinates of 
R 3 x K 3 ) and is equipped with the measure d\(S>d 2 s. Here, d 2 s is the usual rotation invariant measure of S 2 
induced by the Euclidean norm of M 3 and dA is the Lebesgue measure. In the Jaksic-Pillet glued representation 
the corresponding Hilbert space is the antisymmetric Fock space 

The cyclic vector fin is the vacuum of J 7 — ( b 2 ) ■ The representation map ir-jz, of the fermionic C*-algebra V-ji is 
the C*-homomorphism uniquely defined by 

V/ef)i: 7r R (o(/)) = a(g f ) , (6.4) 

where gf € f)2 is given by 



V(p,^MxS 2 : g f (p,#) = \p\(l + e-P p ) 1/2 



f(p&) , P>0, 



f(-p&) , P <0. 



See also (2.16) 



Note that uj-ji is faithful and we thus denote ir-ji (A) and it-ji (V-ji) by A and V-r,, respectively. The weak closure 
of the C*-algebra tt-ji (V-ji) = V-ji is the von Neumann algebra 9JIr := V^. The family {T^} tG R of Bogoliubov 



automorphisms on the algebra V-ji analogously defined as in (2.1 1 uniquely extends to an automorphism group 
of the von Neumann algebra VJI-jz, again denoted by r n := {r^}tea- I n particular, (9Hr,t tc ) defines a W*- 
dynamical system. The thermal state uj-ji is a (/3, r K )-KMS state and so, J7r is cyclic and separating for 9JIr. 
Note that 971 = B(C d ) <g) 9Jt K . 

As above, there is a unique unitary representation of t u whose generator iL-ji satisfies L-jzfl-ji = and 
r^(-) = e ltL (-)e~ ltL for all iel. L K is the second quantization 

L n = dT(p) 

of the multiplication operator by p € R, that is, the operator acting on f) 2 as (pf)(p,'&) — pf(p,i)). 



4. We can now explicitly give the time-dependent C-Liouvillean Ct- By Lemma [2. 1| 

L = L at ® 1^ + ljj, ® L K = (fl- - F at ) ® + lj5. ® dT(p) . 



See also (2.8) 



The time-dependent C-Liouvillean £ t then equals (6.1) with 



W t = r} cos(zat)Hv <S)l n + X^Qi® —= (a(g e ) + a+ (g e )) 



and 



JA^WtA-^J = 77cos( ro t) p at 1/2 g p ^ t /2 ® 1 TC 



A £ f>:t /2 Qipll 2 ® (a(<e-*5/) + « + (%)) 



■ i 

' 4 
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Here, ge(p,{}) := gi(p) (see (2.16)) and dr(lf, 2 ) is the second quantization of lf, 2 , i.e., the particle number 



operator acting on the antisymmetric Fock space Stc := J r _(rj2). 



5. We identify the Hilbert spaces Sat ® {^7?,} C S and Sat and define the orthogonal projector P at E S(S) on 
.flat = Sat ® {Q-r.}- Next introduce the operator C^(t) e B(Sat) defined, for each t G K and any e > 0, by 



(<0 
-ft 



eGo-(L a 



(<0 
at 



ox 



Vi 



(L-e + ie)' 1 



d_ 

OX 



P. 



(<0 



with 



K := Wt - JA 1 / 2 W f A~ 1 ^ 2 J . 



Here, P^t is the Riesz projector related to the eigenvalue e G er(L at ) = <T([£f a t, • ]) of the free C-Liouvillean 
Lot ■ Sat ® {^7?.} -> Sat ® {^-r} defined by 

4t := »^ ® 1«« ~ iJ& 1/2 {H*t <8> l i5K )A- 1 / 2 J . 

Observe that 

L at (AO ) = (Z at A)ft 

for all A G B(C d ) and thus cr(I/ at ) = er(L at ). Note also that £^ does not depend on the time t. 
It is easy to check that the following limit exists 

Ctz := limT^. 

This limit operator C-jz is called the level-shift operator of the atom-reservoir interaction at time i 6 1. See, 
e.g., [271 Eq. (3)]. 

By replacing the terms in £ t containing the coupling constant A by the corresponding level-shift operator 
\ 2 £n we obtain a new generator £[ cff ^ having Sat = Sat ® {^k} as an invariant space. Note further that also 
the subspace 

Sat := {AQq G Sat : ^ G S(C d ), L at (A) = eA, e G {#i - E N ,0, E N - E t }} 

is an invariant space of C-n and Cf for all tel. 

Denote by {U t . s }t> s the two-parameter family of operators on Sat generated by C t . From direct calcu- 
lations 



4 cS \An ) = (zr">*A)n 



(A,»j) 



for all Afi G Sat, where 2[ x ^> i s the adjoint of the Lindbladian £f' v> . Thus, for all A e D 



where {ft*'* ?, '}t>« is the two-parameter family of operators on B(C d ) generated by £ 



Let p{t) := (f£_f )(p at ). Then 



for all A G S(C d ) and thus 



X /flat \ 



We used above that £P^^ = £,f X ' v K Hence, the time-dependent density matrix p(t) satisfies the master equation 

Vt>0: ^p{t)=^ ) p{t), p(0)=p at . 
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6. In order to justify this procedure note that, for all A € 35, 

( P (t),A)^ t = (no, ^(An,,)) , (p at (t), A)^ at = (n , c/ t ( , ' n) (An )) , 

and that in pQ we show, for /o at G 33, 



(off) 7-r(A,r;) 



Co' W ](^o) 



<C ro p|||A| 



for some constant C ro € (0, oo) depending on w but not on the initial state w at of the atom, the observable 
A G 3), and the parameters t, A, and n with |A| sufficiently small. 



6.2 Proof of Theorem H7 



Before starting the proof, we first extend the definitions of A( A ^) kjgh, i) °' 0) kjSj, U k-36), and A( A ^ 



( |4.37[ ) to all eigenvalues e € <r([i/ a t, ■ ]) as follows: 



j^ ) := spanjw^' "' } 



e (J t e+ „, ne{i,...,ni}, n'efi,...,^}} ■ 

m£{-2,-l, 0,1,2} 



The unitary operator U £ is defined by 
whereas 

A (x,n) .- p(o,o) G (A,n) P (o,o) and := ILA^U! 



See, e.g., (3.7), (3.8) and (4.33|. Observe that 



(0,0) _ 2(0,0) _ =(0,0) 



whereas 



Ve G a([ff at , • })\{-w,0,w} : i^ 0) J_ j^ '°>. 



(6.5) 



(A,*7). 



Using these observations we can deduce the spectral structure of the operators A 

Lemma 6.1 (Spectral properties of operators Ai*'^) 
For all e € cr([if at , • ]), 

cr{Ai x ^)\{-ie} C iM-M+ 

and — itn, 0, izu are simple eigenvalues of A^ , A A,,? '' and A^ 7 '', respectively. Moreover, for all e G 
<j([H aU ■ })\{~w,0,m}, 

cr(A^'^) C 

I.e., any eigenvalue of h.^'^ has a strictly negative real part for every e G a([H at , • ])\{ — te, 0, w}. 



Proof. Similar to the proof of Theorem |4.6| (i), we compute that 



fl!:" 



(6.6) 



Indeed, if e G {— vo, 0, nj} then = ^o°'°' ) and the computations for this case are exactly those given ir 



the proof of Theorem 4.6 For e G <r([-ff a t, ■ ])\{ — 0, vo\ the properties (2.10H2.12) of the optical pump 
yield £ p U(o,o) = whereas one can readily check that the Lindbladian Z-n conserves the atomic subspace ^e°'°^. 



These properties lead to (6.6) for e G cr([H at , ■ ])\{— w, 0, nj}. 
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As a consequence, 



a(A^) = a(A^) = -re + <x( A^ | (0 0) ) , 



(6.7) 



where A^ A,?, ^ is seen as an operator acting on fj at , see again proof of Th eorem 4.6 Moreover, he is an 
invariant space of A( A,I? ' g B(f) at ). As explained in the proof of Theorem 4.6 the operator G S(i5 a t) 

is the generator of a relaxing, Markov CP semigroup, see Definition 6.9 In particular, is a non-degenerate 



eigenvalue of A^ X ' V ^ and the corresponding eigenvector is an element of Sj 



(0,0) 
±w 



to 



(0,0) 



As a consequence, from 



(6.7), — ivj, and iw must be non-degenerated eigenvalues of AL , A ,v> and A_£ , respectively. 



Finally, by Theorem 6.10 



Re{a(A( A ^)}\{0} C (-oo,0) 
which, for any e G a([H a t, ■ ])\{— vj, 0, m}, implies that 



Re I a A^ 



= Rela A^ 



\{0}c(-oo,0) 



because of (6.5). 



□ 



The proof of Theorem |4.4| needs further technical results. The next one concerns the stability of the irre- 
ducibility of quantum Markov chains (Assumption [2| under block localization: 

Lemma 6.2 (Stability of Assumption [2] under block localization) 

Let H G B(C d ) be any self-adjoint operator and denote by 1 £ the orthogonal projection onto the eigenspace 
of [H, ■ ] associated with the eigenvalue e G &([H, • ])■ Then, there are rh G N, non-negative real numbers 



{cfl} ■ k i C Rq > rea l num bers {df k } j k ~ t C R, and operators {V^ k '} l k j C B(C d ) such that 



?(f) 



^ l e £nl e = — «[i?Lamb, ■ ] + £<2 j 

6€<7([H,-]) 



where 

and, for any p G .flat, 
Moreover, 



H 



Lamb 



E E E«OT 

:G<r([ff at ,-])\{0} (j,fc)et e £=l 



e e E€H?^r-^^-^r^) 

ee<r([J? at) -])\{0} (j,fe)6i« 1=1 



( U {*$})" = *(<*). 

In particular, £k is fie generator of a relaxing, Markov, CP semigroup satisfying Assumption^ 

Proof. Let {Ej}jL 1 be the eigenvalues of the self-adjoint operator H (M < d) and 

t e := {(i, fc) : Ej - E k = e} C {1, 2, . . . M} x {1, 2, . . . M} . 
For any V G B(C d ) and each eigenvalue 

ee#, ■]) = { J B J --4:i,feG{l,2,...M}} , 
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we define 

V e :=l e (V)= 1[H = Ej] V 1[H = E k ] e 

(j,k)eu 

By construction, note that 

E v. = v- 

ee<r([H,-]) 

In particular, one has 

Fespan{ (J V s \. 

ee<r([H,-]) 

For (i, fc) G t e , e e CT([fl» t , -])\{0}, and e G ct([#, • ]) let 



(6.8) 
(6.9) 



c^ e) := c W 



a j,k — a j,k 



-j,k — "j,k > "i.fc — ™j,k • ( 6 - 10 ) 
By identifying the finite sets {1,2, ...m} x cr([iJ, • ] and {1,2,... to}, we infer from (6.9| and Assumption [2] 
that the family {V^jj. } 3 - & f C £>(C d ) of operators satisfies 



( U 



(6-11) 



{(J,k,i) : cf^jtO} 



For all eigenvectors A ei ,B e2 G /3(C rf ) associated with eigenvalues £i,£2 

€ cr([H, ■ ]) respectively, i.e., l £l (A) = A ei and 1 £2 (B) = £? £2 , let 

D :=A ei B e2 e B{C d ) . 

If D 7^ then we observe that (£1+62) € c([H, ■ ]) and l El + £2 (D) = D, i.e., 13 = D £l+e2 must be an eigenvector 
associated with the eigenvalue £1 + £2. Therefore, by using ( |6.8[ ) and the equality l_ e (V*) = V* , we deduce 
that 



X) l e (Vl e (p)V*) 

e£<r([ff,-]) 



E ^ E Vele(p) £ 
e£«r([H,-]) \eeff([i?,.;i) s£c([H,-]) 

^ ^ v g i e (p)i_ g (n 

g£<r([H,-]) 



for any p G f) a t = /3(C d ). Similarly, 



£ l e (Wl e (p))= ]T V?^p, 

e£o-([H,-]) e£<7([H,-]) 

£ l e (l e (p)V*V) = p y ~* v - 



(6.12) 

(6.13) 
(6.14) 



for any p G f) at . Therefore, by using d3.6k, (3.11)-(3.13), (6.10) and (6.12 1— ( 6.14 ) with the identification of 



the finite sets {1, 2, . . . to} x c r([-ff . ■ 1) and {1, 2, . . . to} we obtain the explicit decomposition of £7?. stated in 
the lemma. By Theorems 



6.7 



and 



6.10 



together with (6.11 1, it is then straightforward to verify that £7?. is the 

□ 



generator of a relaxing, Markov CP semigroup. 

We need more precise information about the behavior of the spectral gap given in Lemma |6.1| with respect 
to the coupling constants 77 and A. This is achieved by using the previous lemma: 

Lemma 6.3 (Behavior of the spectral gap of A^^) 
For all (A, 77) G R x R, 

min{ |Re{p}| : p G a(A^)\{0}} > C ro A 2 
with C ro G (0,oo) being a constant depending on zu but not on A, 77. 
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Proof. We define the function 



g(X, 77) := A" 2 min { |Re {p}\ : p G a(A< A '"))\{0}} 



on the set K\{0} x K. Observe that g(X,rj) only depends on the ratio x := 77/A 2 and is strictly positive, by 
Theorems |6 . 7| and |6 . 10| Indeed, by the proof of Lemma [6. 1| 



g(X, 77) = min \ |Re {p}| : p G a ( ^£ p + £ K ) \{0} 



Furthermore, by Remark 3.2 x£ p /2 + £^ is the generator of a Markov CP semigroup satisfying Assumption 
[2] for any (see also Remark |3.1[ ) and Theorem |6.10| yields that this semigroup must be relaxing. 

By Kato's perturbation theory 26J, for some constants C,c £ (0, 00), g (A,ry) > C whenever 77 < cA 2 , i.e., 

x 1 Y g(x~z, 1) is a strictly positive 



6.10 



when k < c. Using again Kato's perturbation theory :26J and Theorem 1 
continuous function on the interval [c, d] for any finite constant d > c. By compactness of the interval [c, d] , it 
follows that 

min lg{x~%, 1) | x G [c, c']| > . 

So, it remains to prove that g{X,rj) > C whenever 77 > c'A 2 , i.e., when x > d, for some constant C G (0, 00) 
and sufficiently large d < 00. By (6.7) for A, e = 0, note that 

<r(A< 0, ">) C iK . 



Thus, by Kato's perturbation theory 



for the spectrum of x£ p /2 + £7^, the limit 
lim g(x~5, 1) G [0, 00) 



exists and satisfies 



where 



lim ff(x-*,l) > min{|RcM| :pe a(£ TC )\{0}| 



(6.15) 
(6.16) 



J 7? 



e ef<r([i? P ,-]) 



Here, l e denotes the spe ctral projection of [i? p , • ] onto the eigenspace associated with an eigenvalue e G 
o"([flp, ■ ]). Using Lemma 6.2 £7^ is the generator of a relaxing, Markov, CP semigroup satisfying Assumption 
[2] and, by ( |6~16| and Theorem [6~T0 , 

lim s(x~3 5 1) g (0,oo) . 

In other words, for some constants C G (0, 00) and sufficiently large c' > c, <?(x~s , 1) > C for all x > c'. □ 



We now use Lemma 6.3 to obtain norm estimates on the difference of the uniformly bounded semigroups 



{e' 



} Q > and {e 



}* 



>()• 



Lemma 6.4 (Semigroup estimates) 

-For sufficiently small \X\, there are constants C ro , c ro G (0, 00) depending on w but not on X, 77 and a such that 



Va > 



-«A 2 



') 



Proof. Note that the continuous semigroup {e aA< } a >o on fj a t can be represented through the inverse Laplace 
transform of the resolvent of its generator. Indeed, by Lemma 6.1 combined with [5TJ Proof of Corollary 5.15], 

rw+iL 



e« A<A "' , -E= lim \± 



[ + e az f(z - A^)yi --\dz \ 

Jw-iL V Z J J 



.17) 



42 



E 



(z - A (A,?,) ) _1 dz 



(6.18) 



for any w > and with 

— I e az (z-A^)- 1 dz= — 
2iri J 2tti 

\z\=R> \z\=R> 

and R' > sufficiently small. Note that is a simple eigenvalue of A^- 17 ' and the last equality follows from the 
fact that the map 

z^ (1 - e az )(z - hS^)' 1 



is holomorphic near z = 0. By (4.21) and Lemma |6.3| the Riesz projection E associated with the generator 
A( A ' ?? ) is well-defined f or su fficiently small R' <E (0, cX z ) at fixed A ^ 0. Using again the s pectra l properties of 
given in Lemma |6.3[ we can push the integration path of the complex integral in ( |6.17| ) to sufficiently 
small, but a strictly negative real part w' — — cA 2 < as follows: 



lim 



2iri 



w-\-iL 



w—iL 



lim 

L— J-oo 



2m 



c az (^(z-A^)- 1 - dzj 

e az y(z - A^'^) _1 - dz 



+ 



2ni 



-cX 2 +iL 

-c\ 2 -iL 



lim 

L— >oo 



\z\=R> 
I 2m 



1 

2m' J 



dz 



-cA +iL ctz 



-c\ 2 ~iL z 



( z _ A (A,r,))-l A (A,r,) dz L +E_ 1 . 



Note that c £ (0, oo) is a sufficiently small constant depending on vo but not on on A, rj, and a. By (6.17), it 
follows that 



lim <^ A" 



2iri 



L ~a(ix—c)\ 2 

l (ix - c) 



{{ix - c) - A" 2 A( A '"))- 1 A^da; I (6.19) 



for all a € Kg". Observe that we ha ve a d ditio nally used an obvious change of variable in the last equation to 
extract t he fac tor A~ 2 . Indeed, by f(^-(|6j} and Assumption [j] (| 77 1 < CA 2 ), we have ||A( A ' r ')|| = C(A 2 ) and 
the limit (6.19) yields the upper bound 



E|| 



<e~ caA lim 

L— too 



2m 



ix — c 



{{ix-c)- A^A^)" 1 dx\ (6.20) 



for all a £ Kq . Hence we need to bound the integrand of the last integral by some integrable function not 
depending on the parameters A, r\. To this end, first observe that 



A -2 A (A,„) = u* ( Cn + \ Uq 



.21) 



with the ratio x := ry/ A 2 € [— xo,>fo] for some fixed kq £ (0, 00). See (4.37), Theorem 4.6 (i) and Assumption 
[T] Therefore, we define the map g :R x [— xq, xq] — > K by 



< 00 . 



.22) 



g{x,x):= {{ix - c) - £ n - x£ p ) 1 
This function is well-defined for sufficiently small c > because of Lemma |6.3| Moreover, since 
\g{xx, kx) - g(x 2 , x 2 )| 

< |3(a; 2 ,x 2 )|(|| a;i -x 2 || + !|x 1 -x 2 ||||£ p ||) , 



43 



-(n+l) 



the function g is (locally Lipschitz) continuous on its domain of definition. By compactness, for any L £ Rq~ 
there is a constant C £ (0,oo) depending on w but not on A, rj such that 

sup g(x, h) < C . (6.23) 

(x,x)E[ — L,L] X [-3*0,Xt)] 

On the other hand, the Neumann series 

oo 

((ix -o)—Cn — x£ p ) _1 = y^^jix - 
implies that, for sufficiently large x » 1 and h £ [— >c , *e ], 

oo 

~ 9 (x, *)<J2 ii« - cir (n+1) {ii^ii + ii£ P iir 

n=0 

because £7^ and £ p are bounded operators on f) at . As a consequence, there is a sufficiently large constant 
C £ (0, 00) depending on w but not on A, r\ such that 



Vi e R, x £ [— xq, xq] : <?(a;, >f) < min < C, 



isc — c 



(6.24) 



By (6.20H6.24), we conclude the existence of a constant C £ (0, 00) depending on w but not on A, r\ such that 

11 qA< a <»> 



Va > : 



Ell < Ce 



-caX- 1 



(6.25) 



In the same way we obtain (6.20), we derive the upper bound 

caX2 lim 



< e" 



Here, 



E 



L— voo 

1 2tti 



2-7TZ 



f 



— L 11^ " c ll 



((ix-^-X^P^G^P^)- 1 do; I. (6.26) 



e a, (z _ Po (A^) G (A,r,)p o (V)))-l dz 



.27) 



is an operator associated with the generator P^'^ G^ X ' V ^ Pq X ' v \ The latter is well-defined for sufficiently small 
R' £ (0, cA 2 ) at fixed A ^ 0. 

Namely, for |A| <C 1 and some constant C £ (0, 00) depending on w but not on A, 77, observe that 

||P o (V7) G (V7)|| <CA 2 and ||G (a, ''^Pq '°'|| < CA 2 , 



using P^ ' )^ ' ) = 0. Hence, combining these last upper bounds with (4.30) and the triangle inequality we get 

= ||P (A ^G (A '" ) P (A ' I ' ) - P (0 ' 0) G (A ^P ^ 0) || < CA 4 (6.28) 

wit h C £ (0, 00) n ot de pending on A, r\ and ru. In particular, the spectral properties of A' A ' r " given in Lemma 
together with (6.28) implies that E is well-defined for a sufficiently small, but strictly positive R' £ (0, cA ) 

E - E|| - f A 1 -i \\{z - A^)- 1 !!!!^ - P^G^P^y^ldz < CA 2 (6.29) 



6.3 



and satisfy 



\z\=R> 



for some constant C £ (0, 00) not depending on the (sufficiently small coupling) constants |A|,|r?| 
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To 



prove the inequalities 



(6.26) 



and 



(6.291, 



note 



that 



the 



(non- 



degenerate) eigenvalue is the unique element of the spectrum of A( A ' r ') within the disc of radius R' , 
by Lemma 6.3 By (6.28) and Kato's perturbation theory [26 of discrete eigenvalues, there is a unique 
(non-degenerate) eigenvalue of p^^ q( x ^) p^' 71 ^ anc [ no other element of its spectrum within the disc of 
radius R', provided that |A| is sufficiently small. In fact, this eigenvalue is also because G^^ is a closed 
operator and the constant function Slo € .flevo satisfies G^ x,r, ^*Qo = for the adjoint G^'^* of G^' 17 ' (implying 
that is also an eigenvalue of G (X ' v) and thus of P (A,, ' ) G (A '' ?) P (A ^' ) ). Hence, the map 

z ^ (i _ e «)( z _ p(^) G (\n) P ^v)yi 



is holomorphic near z = 0. Consequently, 

1 



E = 



2iri 



(z-P^G^P^^dz 



\z\=R> 



is the Riesz projection onto the eigenspace associated with the eigenvalue of p( x,ri ) q(\,v) p^'ri an( ^ j-,y (6.18), 
the first inequality of (6.29) follows. 



Now, using ( 6.28 ) , the fact that g(x, x) is uniformly bounded, and the second resolvent equation for sufficiently 
small |A|, we get the upper bound 

c ) _ X^P^G^P^)- 1 < C g(x, v /X 2 ) 



for some constant G € (0, oo) depending on w but not on A, r\. Therefore, similar to (6.25), we infer from (6.26) 
and properties of the function g that there is a constant G S (0, oo) depending on w but not on A, r\ such that 



Va > 



\e aP o 



Ell < Ge 



-caX- 1 



Combining this with (6.25) and (6.29), one obtains the statement of the lemma. 



□ 



We now conclude by the proof of Theorem 4.4 that is, we show, for sufficiently small |A| and any s € (0, 1), 
that 

lexp (aA^>) - exp (aP^' n) G^ P^ v) ) II < G ro , e |A| 2(1 - £) , (6.30) 
where the constant G roj£ G (0,oo) depends on w,e but not on A, ij, and a. 

Proof of Theorem 



u 



Note 



that 



the 



semigroups 



and 



{e aP o ' 7) G (A ' r,) P ( ,!> ]. Q > are uniformly bounded in A, 77 as the first one is a CP semigroup and the sec- 



ond one is the restriction of {T a } a >o, see (4.7). Thus, DuhamePs formula yields the inequality 



p " || < aCWP^G^P™ - A( A '")|| 



for some constant G € (0, 00) not depending on A, 77, 137, and a. By (6.28), it follows that 



eO pt» fl ) 0Mp M ||<cf * 



.31) 



(6.32) 



with G ro G (G, 00) depending on w but not on A, rj, and a. Finally, we infer from (6.32) and Lemma 6.4 that 



„aA< A -"> 



||e~" -e 
< G ro min {aA 4 , A 2 + e~ ac ^ A * } 



(6.33) 



< |A| 2C1 ^ e) (1 + |A| 2£ + |A|- 2(1 "^ } < 



,|A|" 



for any e € (0, 1) and we obtain (6.30) for sufficiently small |A|. 



□ 
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6.3 Completely positive (CP) semigroups 

In the theory of open quantum systems, one is usually interested on the restricted dynamics of some small 
quantum object interacting with macroscopic reservoirs. This restricted time-evolution is described by a map 
C on the set of density matrices of the small system. See, for instance, [21 Section 1.2.1]. Properties of such 
maps (cf. [31 Section 1.2.2]) motivate the definition of the class of completely positive (CP) operators: 

Definition 6.5 (Completely positive (CP) maps) 

A positive map C G B (B (X)) acting on the set B{X) of bounded operators on a Hilbert space X is called 
completely positive (CP) if the extended map C ® ls(c™) remains positive for all n G N. If C (lx) = ^-X> then 
the operator C is called a unital map. 

Completely positive (CP) semigroups are defined as being the semigroups which are CP maps for all times: 
Definition 6.6 (Completely positive (CP) semigroups) 

A semigroup {Ct}t>o C B (B (X)), with X being a Hilbert space, is CP if the map C± is CP for all t G Rq~. If Ct 
is unital for any t G Rj, then we call {Ct}t>o unital. 



From now on and until the end of Section 6.3 X is always a 
n-dimensional Hilbert space. We denote by B2 (X) = B (X) the Hilbert space of Hilbert -Schmidt oper- 
ators with scalar product 

(A,B) B2(x) :=Tr x (A*B), A,BeB 2 (X). 

In the special case where a semigroup {C t }t>o C B (B 2 acts on B2 (X), we can define its (unique) adjoint 
semigroup {C t *}t>o C B {B2 (X)) as usual via the equations 

Vi>0: (C:(A),B) B2{x) = (A,Ct(B)) B2 ( X ), A,BeB 2 (X). 

Note that if the CP semigroup {Ct }t>o is unital then {C t *}t>o is CP and preserves the trace. A CP semigroup 
{Ct}t>o is called Markov CP semigroup if it preserves the trace. Generators of Markov CP semigroups {C t }t>o 
and of their adjoint groups {C t *}t>o can be characterized in the finite dimensional case as follows (cf. [42"1 
Theorem 2.2] and gl Theorem 2]): 

Theorem 6.7 (Generators of Markov CP semigroups I) 

Let dim A" = n G N. The operator L G B (B2 (X)) is the generator of a Markov CP semigroup {C t } t >o if and 
only if 

L (p) = -i [h, p] + - £ { [V, , pV*} + [V jP , V*] } , AeB 2 (X) , 

3 

where h = h* G B2 (X) and {Vj} C B (X) is a family of operators with ^2j~Vj~Vj G B{X). Additionally, the 
adjoint semigroup {C t *}t>o is in this case the unital CP semigroup with generator 

L* (A) = i [h, A] + 1 £ { [V* , A] V, + V* [A, V,]} , A e B 2 (X) . 

3 

Generators of CP semigroups are also called Lindbladian or Lindblad (-Kossakowski) generators. 
A more compact characterization of generators of Markov CP semigroups is given by [571 Sect. 4.3]: 

Theorem 6.8 (Generators of Markov CP semigroups II) 

Let dim A" = n G N. The operator L G B (B2 (X)) is the generator of a CP semigroup {Ct}t>o if and only if 
there is a completely positive map S G B (B (X)) and an operator A G B (X) such that 

L = + ^ + S . 

Such a CP semigroup is Markov if and only ifL* (1-x) = 0. 



4G 



The relaxing property of Markov CP semigroups {Ct} t>0 C B (Bi (X)), which is crucial for our analysis, is 
defined as follows: 

Definition 6.9 (Relaxing semigroups) 

A Markov CP semigroup {Ct}t>o C B (B2 (X)) is called relaxing if there is a unique trace-one positive p x G 
£>2 (X), i.e., a density matrix p^, such that, for any density matrix p € B2 (X), 

lim C t (p) = Poo . 

t— >oo 

In other words, a relaxing, Markov CP semigroup has a unique invariant equilibrium state. Moreover, this state 
can be approximated by the density matrix Ct (p) for large times and any initial state with density matrix p. 
Spohn [301 Theorem 2] gave in 1977 a characterization of relaxing semigroups which turns out to be very useful 
in our context: 



Theorem 6.10 (Condition for a Markov CP semigroup to be relaxing) 

Let dim A" =n£N. Let {Ct}t>a C B {B% (X)) be a Markov CP and Cq semigroup with generator L given by 
Theorem 6.1 If the space spanned by the family {Vj} satisfies 

span{Vj} = span {V*} C B 2 (X) 

and the bicommutant 

{Y 3 }" = B(X) , (6.34) 
then {Ct}t>o is relaxing. In particular, is a non-degenerated eigenvalue o/L and 

max {Re {w} | w € a (L) \{0}} < 0. 

As explained after Assumption [2j the condition (6.34) is a non-commutative version of the irreducibility of 
classical Markov chains. 



6.4 The operator Q3 p .rc and the conservation of positivity 

1. We first give an elementary example of an operator Q3 Pi -r, defined by (4.44) which does not preserve positivity. 
Let d — 3, N = 2, zu = 1 and 






in the canonical orthon orma l basis of C 3 . Choose the coupling functions it such that gi(e) > for 
and e e {-1,0,1}, see \2.m\. Then, the family {V$}j,k,l C B(C 3 ) of operators defined by (3.9 1 for t E {1,2} 



and j, k € {1, AT} satisfy Assumption 



defined by (4.33) are eigenvectors of 



i 2 With this choice, the operators W^jl' 
Cti with eigenvalues 7 1 ,7 2 ,7i,72 G C\{0}, 



e{i,2} 



W (N,2) W (TJT TJ/Cl.l) c fi 
' ^(1,1) ' W {N,iy W (N,2) G ^ at 

respectively. By setting 




we infer from (5.2) that 























-( 




v 





J 





in the canonical orthonormal basis of 1 
10 





2Re{ 7 r 1 } 

7T" 1 +7 2 _1 

Indeed, note that 





7r 1 +7 2 ~ 1 

2Rc{ 72 - 1 } 



ftp I 




= w, 



(1,1) 



w, 



(N,2) 

(1,1) ' 




K = W, 



(1,1) 



w, 



(1,1) 

(JV,2) 
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and 
Thus, 



VV (1,1) P ~~ Wl) + ^(JV,2)' "P^(JV,j) — K Vj) (JV,J) ' J l- 1 '^!- 



B- 




iv.i o o o =-( 7r 1 ^ + 7 -<^)/ i ; 



Let 



- ^fl*® + <;J) - 72 1 (^iv,'i) + <j?) 



D-deti ^^r 1 } T^+Ti" 1 , 

^- det l TI^+T^ 1 SRc^l ) G 



be the determinant of the lower 2x2 block diagonal part of the above self-adjoint matrix. One can clearly find 
coupling functions for t E {1,2} such that |gi(l)| 2 = |g2 (1) | 2 , but 

VV( gl )(l)^VV(g 2 )(l) . 

This implies that 

Re^ 1 } - Rc^- 1 } and Im^ 1 } ^ Im{ 7 ^ 1 } . 

However, in this case 

^ = -(lm{ 7 r 1 }-Im{ 72 - 1 }) 2 <0 

showing that Bjy.i is not a map from B + (Hi) to B + {Hn) and hence 23 p ,-r is in general not the generator of a 
semigroup on D which preserves positivity. 

2. Note that under certain conditions the weighted stimulated transition rates defined via the operator 

/>oo 

Sp.K := / c- s£K £ p e s£K £ P ^D ds 
Jo 

/>oo 

= / e- 8fl *£pf#e' fi *£pP»cfceiB(i5 Bt ) 
Jo 

are well-defined. This is the case, for instance, for strong decoherence, i.e., if 

min(Re{cr(£ K P s ,)}) > max(Re{cr(£ K P^)}) 

as in this situation 

|| e - S £ K£pP ±|| Ue^fipP^II < Ce~ cs 

for constants C, c e (0, oo) and all t > 0. Observe that the norm ||e _ * £TC || exponentially grows as t — > oo, 
generally. In this case we obtain, moreover, that 

In particular, 2$ Pl rc is the generator of a positivity preserving semi-group on 55. 

Therefore we could wonder whether or not 25 p ,-r also generates a positive semigroup if the system shows 
strong decoherence. Unfortunately, this turns out to be not the case. This can be seen in a simple variation of 
the example above: By introducing a third atom-reservoir interaction term in the model above with coupling 
satisfying | g3 (0) | 2 > and 

q 3 
Qs ~ | | , q 3 G 




48 



the diagonal part £-r,Po is not changed whereas ^-rP^ is shifted by 

— I §3 (0)1 2 -F® • Hence the strong decoherence regime can be always attained 

for sufficiently large q 3 . As the eigenvalues 7i)72>7ii72 OI ^(^lj^^^fi) ' 

^Wil'^Wj) are simultaneously shifted by — |g 3 (0)| 2 , P remains negative for all g 3 G R. 

3. Nevertheless, Assumption [3] can be verified for a large class of models. For instance, this condition is always 
satisfied if the 1st and the Nth. energy levels of the atomic part of the model are non-degenerated. This is a 
special case of the following theorem: 

Theorem 6.11 (23 p ,-r, as transition rates) 

Assume that Pn,i&r.Pn,i = £jv,i-Pjv,i for some £ N1 G C. Then Assumption^ holds. 
Proof. Using (3.6 1 and (3.11 )-( 3.13[ ), we observe that 

InaPi,n (p) = (t Ntl PN,i (p))* = {Pn,i£kPn,i (p)T = Pi,n&rPi,n (p) 

for any self-adjoint p = p* G -f) at . Since any element of $j a t is a sum of self-adjoint elements and £>tz as well as 
Pi, at and P\ t M are linear operators, one obtains that 

P\,N&rP\,N = £,N,lPl,N ■ 

Using the notation £ := £ N i, we get that 

Pv,l (^lran(P^)) P N,1 = £ 1 Pn,1, Pl,N (^lran(F^)) Pl,N = I 1 Pl,N ■ 



Thus, by \5.2\ , 

Bat,! = — £ 1 Pn ,n1t* p Pn ,ih\yPi,i — £ PN,Nh^Px,N^ v Pi,i 



Using the properties of h p and ^4 £? = ^,4, 



Similarly, one also finds 



On the other hand, we have 



_ 2Re(£) 
B JV.i |^p— ^ £p • 



_ 2Re(£) 
B i.iv |£p~ £p 



Biv,jv — 6 hj>h^>PN,N + 6 1 PN,Nhph* p , 
B M = r'^Pi.i + r'Pi.i^Jp • 
Therefore, by (5.3), we can split the operator Q3 p k on the subspace I? as follows: 

®p,k = ^ + ^!+S (6.35) 

with 

a := (« + kk) + + w 

- r 1 (m; + Kh p ) g 6(c d ) 

and 

2 := B 1:A r + Bat,! G B(H at ) ■ 
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Since Zn is the generator of a relaxing Markov CP semigroup on Sj at (Remark 3.2), Re(£) < and S is thus a 
completely positive map. Clearly, 



(^+^ + ~)* (i) = (a^+^+s) (i) = o 



Therefore, we infer from ( |6.35 ) 
Markov CP semigroup on f) at . 



and 
As 



6.8 



that 25p,7?. generates a 
an invariant subspace of this 
generator, S Pi k generates a Markov semigroup on T> which preserves positivity. □ 



Theorem 
V is 
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